Brane-induced Skyrmion on S^: baryonic matter in holographic QCD 
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We study baryonic matter in holographic QCD with D4/D8/D8 multi-D brane system in type 
IIA superstring theory. The baryon is described as the "brane-induced Skyrmion", which is a 
topologically non-trivial chiral soliton in the four-dimensional meson effective action induced by 
holographic QCD. We employ the "truncated-resonance model" approach for the baryon analysis, 
including pion and p meson fields below the ultraviolet cutoff scale Mkk ~ IGeV, to keep the 
holographic duality with QCD. We describe the baryonic matter in large Nc as single brane-induced 
Skyrmion on the three-dimensional closed manifold with finite radius R. The interactions between 
baryons are simulated by the curvature of the closed manifold S^, and the decrease of the size of 
represents the increase of the total baryon-number density in the medium in this modeling. We 
investigate the energy density, the field configuration, the mass and the root-mean-square radius of 
single baryon on as the function of its radius R. We find a new picture of "pion dominance" near 
the critical density in the baryonic matter, where all the (axial) vector meson fields disappear and 
only the pion field survive. We also find the "swelling" phenomena of the baryons as the precursor 
of the deconfinement, and propose the mechanism of the swelling in general context of QCD. The 
properties of the deconfinement and the chiral symmetry restoration in the baryonic matter are 
examined by taking the proper order parameters. We also compare our truncated-resonance model 
with another "instanton" description of the baryon in holographic QCD, considering the role of 
cutoff scale AIkk- 
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I. INTRODUCTION 

The concept of "holography" is first introduced by 
D. Gabor in 1948 [1] as a new technique of the optical 
physics to playback the three-dimensional information 
onto the two-dimensional plate as a hologram. In 1997, 
this concept of holography comes to get a new appear- 
ance in the framework of the superstring theory as the 
duality between two theories belonging to the different 
spatial dimensions. It is first proposed by Maldacena [2] 
as the AdS/CFT correspondence between AdS^ x su- 
pergravity and A/'=4 SUSY Yang-Mills theory through 
D3 brane in type IIB superstring theory. In more gen- 
eral point of view, an essential element of the holography 
is "Dp brane" as the (p-l-l)-dimensional membrane in the 
ten-dimensional space-time. In fact. Dp brane appears as 
the soliton, i.e., the condensed object of the fundamen- 
tal strings. The Dp brane has two important aspects as 
follows: {p + l)-dimensional gauge theory appears on a 
surface of the Dp brane, and {{p + 1) + l}-dimensional 
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supergravity appears around the Dp brane. (The italic 
"i" denotes the radial dimension with non-trivial cur- 
vature around the Dp brane, indicating the existence of 
the gravity). Actually, the concept of holography indi- 
cates the duality between the [p -\- l)-dimensional gauge 
theory without the gravity and {{p-\- 1)-^- l}-dimensional 
supergravity mediated by the Dp brane, and the gauge in- 
teraction as a "hologram" on the surface of the Dp brane 
is to give the supergravity as a "vision" in the extra- 
dimension. 

One of the most essential properties of the holography 
is the "strong-weak duality" between the gauge theory 
and the supergravity: the coupling strengths are trans- 
versely related with each other. Therefore the holography 
provides a remarkable possibility that non-perturbative 
aspects of one side can be analyzed by the other dual 
side just with the tree-level calculations. Then, if we find 
the special configurations of D branes reflecting QCD on 
their surfaces, non-perturbative aspects of QCD can be 
successfully examined from the tree-level dual supergrav- 
ity side. This is the strategy of the holographic QCD. 

There exist several trials to find the special configu- 
rations of D branes reflecting QCD. Eventually, in 2005, 
Sakai and Sugimoto succeeded in constructing QCD with 
massless quarks and gluons from the fluctuation modes of 
the open strings on the D4/D8/D8 multi-D brane config- 
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urations in type IIA superstring theory [3] , called Sakai- 
Sugimoto model, which is one of the most realistic model 
of holographic QCD. By using this model, many phe- 
nomenological properties of mesons belonging to the non- 
perturbative aspects of QCD like meson mass spectra, 
hidden local symmetry [4], vector meson dominance [5], 
KSRF relation [6], GSW model [7], etc., are successfully 
derived from the tree-level dual supergravity calculations. 
In this sense, holographic QCD is often regarded as the 
"unified meson theory" . On the other hand, baryon de- 
scription is not straightforward in this approach since the 
classical supergravity is found to be dual with the strong- 
coupling "large- iVc" QCD, where baryons do not directly 
appear as dynamical degrees of freedom [8]. 

In our previous work, we gave the first study of the 
baryon as a non-trivial topological soliton in the four- 
dimensional meson effective action derived from holo- 
graphic QCD. We call this topological soliton as a "brane- 
induced Skyrmion" [9, 10]. Especially we included pions 
and p mesons appearing below the Kaluza-Klein mass 
scale Mkk ^ IGeV, which is often called the "tnmcated- 
resonance model" for the baryon analysis. 

Actually, Mkk plays the role as the ultraviolet cutoff 
scale of the holographic approach. In fact, there appear 
infinite number of "non-QCD modes" with mass scale 
~ O(Mkk) in holographic QCD like gluinos and Kaluza- 
Klein modes. In this sense, the duality with QCD could 
be maintained below Mkk as the ultraviolet cutoff. 

The appearance of certain cutoff scale Mkk should be 
essential for the holographic approach to be dual of real- 
istic QCD with confinement and chiral symmetry break- 
ing as the non-SUSY natures. In the holographic model 
with D4/D8/D8 multi-D brane system, D4 branes are 
^i-compactifiled with Mkk scale, to give the complete 
SUSY breaking and its resulting non-conformal natures 
of QCD like finite string tension and chiral condensate. 
These considerations suggest that Mkk should be so 
much respected, giving our truncation of meson reso- 
nances at Mkk for baryon analysis (See, Sec. II for de- 
tails). 

Recently a baryon is also described as an "instan- 
ton" on the five-dimensional gauge theory of D8 branes 
with D4 supergravity background [11-15]. The instan- 
ton is introduced before the mode expansion of the five- 
dimensional gauge field into mesons, so that the baryon 
as the instanton is to be composed by the infinite num- 
ber of color-singlet modes with mesonic quantum num- 
ber even above Mkk scale. However, we consider that 
such color-singlet modes above Mkk might not directly 
correspond to physical mesons in QCD, because the du- 
ality with QCD is mentained below Mkk- Furthermore, 
there also exist infinite number of other non-QCD modes 
above Mkk, which could also affect the baryon proper- 
ties if they were included. Therefore, in contrast to in- 
stanton models, we severely respect the cutoff scale Mkk 
and truncate the meson resonances at Mkk, to keep the 
duality with QCD. More comprehensive discussions with 
instantons are summarized in Sec. VII . 



In present work, we newly consider the extension of the 
holographic model to dense QCD. Due to the non-abelian 
nature of QCD, various realizations are expected in the 
vacuum itself with finite temperature and density, called 
"QCD phase diagram". Up to now, interesting phase 
structures are proposed by using some low-energy effec- 
tive theories of QCD, e.g., confined phase with mesons 
and baryons, deconfined phase with quark-gluon plasma 
(QGP), chiral symmetry broken phase with mass gen- 
eration [16], color superconductivity as diquark conden- 
sation [17-21], etc. In fact, some wisdom about QGP 
gives the insight for the early universe just after the Big 
Bang [22]. Furthermore, the possible QCD phase tran- 
sitions in the core region of neutron stars could affect a 
lot of their macroscopic features like moment of inertia, 
angle velocity, and breaking index [23]. There also ex- 
ist several experimental projects to search the QGP in 
the ultra-relativistic heavy ion collisions in RHIC (Rel- 
ativistic Heavy Ion Collider) at BNL, and LHC (Large 
Hadron Collider) at CERN. There will also appear rela- 
tively low-energy collision experiments to make the low- 
temperature high-density object in FAIR (Facility for 
Antiproton and Ion Research) at GSI, giving some knowl- 
edge about the core region of the compact stars. With 
these background, it should be urgently important to 
make clear the structure of QCD phase diagram more 
explicitly from QCD itself with the rich help of experi- 
mental data, which will eventually bring about the fun- 
damental understanding of our whole nature. 

There exists the lattice QCD numerical study as the 
first principle calculation of the strong interaction. How- 
ever, because of the "sign problem", its applicability 
is severely restricted near the zero-density at finite- 
temperature regime of the wide QCD phase diagram (For 
some review, see, Ref. [24]). Therefore, if one succeeds 
in the extension of holographic approach to the dense 
regime, it should give a new analytical tool for non- 
perturbative aspects of the finite density QCD, where 
the holography provides the duality between the strong- 
coupling gauge theory and the weak-coupling supergrav- 
ity. This is the main aim of our study. 

In this work, we consider the baryonic matter in holo- 
graphic QCD as the extension of the holographic ap- 
proach to dense QCD. Especially, we treat the baryonic 
matter with large-A"c because the holographic QCD is 
derived as a large-iVc effective theory. As the general 
property of large- iVc QCD [8], the static baryon mass 
is proportional to 0{Nc), so that its kinetic energy be- 
comes 0{N~^). There also exist the quantum effects 
like the zero point quantum fluctuation energy Eq and 
also the baryon mass splitting Am within the baryonic 
matter, while these correspond to the higher order con- 
tributions of the 1/A^c expansions as Eq ~ 0{N^) and 
Am ~ 0{N~^) [25]. Such large- A^c countings indicate 
that, for sufficiently large Nc, the kinetic energy and 
the quantum effects within the baryonic matter can be 
suppressed relative to the static mass, and the baryonic 
matter comes into the static Skyrme matter. Such static 
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Skyrme matter was first analyzed by Klebanov [26], by 
placing the Skyrme soliton solutions periodically along 
the three-dimensional cubic lattice, which would corre- 
spond to the nuclear crystal in the deeper interior of 
neutron stars. 

Since the cubic lattice treatment is rather cumbersome, 
we employ a mathematical trick to analyze such static 
Skyrme matter, proposed by Manton and Ruback [27]. In 
order to represent some high-density state of the multi- 
Skyrmion system on the three-dimensional flat space R'^, 
single Skyrmion is alternately placed on a surface of the 
three-dimensional closed manifold with a finite radius. 
In fact, the multi-Skyrmion system on R'^ and the single 
Skyrmion on can be related with each other through 
the compactification of the boundary of a unit cell on 
R'^ shared by one Skyrmion. The interactions between 
baryons on R^ are simulated by the curvature of the man- 
ifold S^, and decreasing the radius of represents the 
increase of the baryon-number density in the medium in 
this modeling. Actually, by taking such mathematical 
simplification, one could avoid some complicated analy- 
sis like Monte Carlo simulations on the three-dimensional 
cubic lattice [26] , and he can get some physical intuitions 
as for the baryonic matter qualitatively and even quan- 
titatively [27]. Therefore, by placing the single branc- 
induced Skyrmion on the closed manifold S^, we try to 
analyze the typical features of baryonic matters in holo- 
graphic QCD. Especially, in this analysis, the roles of p 
mesons in the dense baryonic matter will be examined in 
detail from the holographic point of view. 

There exist a lot of works about the extension of 
the holographic approaches to dense SUSY QCD (See, 
Ref. [28] and references therein). With the bottom-up 
construction of the AdS/QCD models [29-36], it was ap- 
plied to dense QCD by introducing the bi-nuclcon con- 
densate [37]. After the discovery of the Sakai-Sugimoto 
model [3] as one of the most reliable holographic top- 
down approaches to non-SUSY QCD with masslcss quark 
flavors, there exist several proposals about the exten- 
sion of this model to dense QCD. For example, in 
Refs. [38, 39], the baryon chemical potential /xb is intro- 
duced by the asymptotic value of a U(l) gauge field on 
the D8-D8 branes as Vq = —ins/Nc, similarly to the in- 
troduction of the chemical potentials to the chiral pertur- 
bation theory by promoting the global chiral symmetry 
to local gauge one [40, 41]. In these holographic analy- 
ses, the density dependence of the "local" properties of 
mesons and baryons like their masses and coupling con- 
stants, and also the phase structure of QCD have been 
successfully discussed. Now, in our paper, we treat the 
baryon as a "non-local" solitonic object as the Skyrmion 
in the Sakai-Sugimoto model, and we analyze the dense 
QCD by the Skyrme matter. Therefore, adding to the in- 
formations about the local natures of hadrons and QCD 
phase structure, we can extensively examine the internal 
structure of the baryon like the size and the field config- 
urations. 

Here we show the organization of this paper and its 



brief summary. In Sec. II, we overview the holographic 
derivation of the four-dimensional meson effective action 
with the pion and p meson fields, reemphasizing the role 
of cutoff scale Mkk in the holographic framework to be 
dual of QCD. In Sec. Ill, we analyze the properties of 
baryons and baryonic matter in holographic QCD. In 
Sec. Ill A, wc introduce the concept of the branc-induced 
Skyrmion on the three-dimensional flat space R"^ . Then, 
in Sec. IIIB, we describes the baryonic matter as the 
system of single brane-induced Skyrmion on the three- 
dimensional closed manifold S^. We derive the expres- 
sion of the hedgehog mass and Euler-Lagrange equations 
for the pion and p meson fields as the brane-induced 
Skyrmion on S^. Sec. IV is devoted to the numeri- 
cal results and their physical interpretations about the 
baryon nature in dense QCD. In Sec. IV A, the baryon- 
number density dependence of the energy density and 
field configuration profiles of single baryon are discussed. 
In Sec. IV B, we propose a new striking picture of the 
"pion dominance" near the critical density, i.e., all the 
(axial) vector meson fields disappear and only pion field 
survives. In Sec. JVC, the baryon-number density de- 
pendence of the mass and root-mean-square mass radius 
of single baryon are analyzed. Wc find some non-linear 
increase in the size of the baryon near the critical den- 
sity as a "swelling" phenomena. In Sec. IV D, we ex- 
plain the swelling mechanism in the general context of 
QCD, and consider its effects on the stability of A^-A 
mixed matter. In Sec. V, we examine the features of 
the dclocalization phase transitions and the chiral sym- 
metry restoration by choosing proper order parameters, 
through which the relations between deconfinement and 
chiral symmetry restoration are reconsidered. In Sec. VI, 
we calculate the critical densities of the phase transitions 
in the physical units with the experimental inputs for the 
pion decay constant f-Tri— 92.4MeV) and p meson mass 
mp{~ 776.0MeV). Wc find the critical density pB — 7po 
in the holographic approach. Through all of the sections, 
by comparing the branc-induced Skyrmion and standard 
Skyrmion without p meson fields, the roles of the vector 
mesons in the dense baryonic matter are examined from 
the holographic point of view. Sec. VII is devoted to sum- 
mary and outlook. In this final section, we compare our 
truncated-resonance approach with another "instanton" 
description of baryons [11-15] in the holographic QCD, 
paying attention on the role of cutoff scale Mkk- 



II. MESON EFFECTIVE THEORY FROM 
HOLOGRAPHIC QCD 

In this section, we overview the derivation of the 
four-dimensional meson effective action from holographic 
QCD with D4/D8/D8 multi-D brane system in type IIA 
superstring theory, called the Sakai-Sugimoto model [3]. 
The meson effective action is derived without small am- 
plitude expansion to discuss a baryon as a large ampli- 
tude chiral soliton. (For more comprehensive derivations. 
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see our previous paper [9]). Here we especially emphasize 
the roles of cutoff scale Mkk in holographic QCD. 

First, we review the construction of D4/D8/D8 multi- 
D brane system and also the supergravity description of 
D4 branes as the holographic dual of QCD. As the first 
step, Nc sheets of D4 branes are prepared to construct the 
gluon sector of QCD. The D4 branes are S'^-compactified 
along one extra-dimension with a radius as the inverse of 
the Kaluza-Klein mass scale Mkk- There appear ten in- 
dependent fluctuation modes from the open strings on 
the surface of D4 branes, i.e., gauge fields -4^=o~3, scalar 
fields A4 and $^=5^9, and also their superpartners as 
fermions. By imposing the anti-periodic boundary con- 
ditions for all the fermions along the 5'^-compactified di- 
rection, they acquire large masses 0{Mkk)- Then 
supersymmetry (SUSY) is completely broken and, due 
to the radiative corrections, all the scalar fields A4 and 
$i also get large masses ^ 0(A'/kk)- Because of the 5^- 
compactification, there also appear the infinite number of 
the Kaluza-Klein modes with masses ~ O(Mkk)- There- 
fore, below the Mkk scale, only massless gauge fileds 
appears. In this sense, the system of Nc sheets of D4 
branes with the S'^-compactification can be viewed as the 
U(iVc) Yang-Mills theory below Mkk scale, correspond- 
ing to the pure gauge sector of QCD. As the next step, 
Nf sheets of D8 and D8 branes are added to introduce 
the massless quark flavors of QCD. D8 has opposite chi- 
rality relative to D8, providing U(-/V/)l x V{Nf)j^ chiral 
symmetry in this model. From the fluctuation modes 
of open strings between D4 and D8 (D8), there appear 
massless chiral fermions as quarks in QCD. As a whole, 
massless QCD appears as the "hologram" on the surface 
of D4/D8/D8 branes. 

Then we shift into the gravitational description of D 
branes from the extra-dimensions. The D brane is origi- 
nally introduced as the fixed edges of open strings with 
Dirichlet boundaries. This also indicates that, through 
the "open-closed duality" for the fundamental strings, 
the D brane can also be regarded as the source of closed 
strings, giving the graviton in the extra-dimensions out- 
side of the D brane. In this sense, the D brane can be 
identified as a highly gravitational system, i.e., the "black 
brane" , allowing the gravitational description from the 
extra-dimensions. In fact, the mass of the D branes is 
proportional to its sheets number, so that, by assuming 
Nc ^ Nf, only D4 branes can be represented by the grav- 
itational background and D8 (D8) branes are introduced 
as the probes called "probe approximation" , which cor- 
responds to the quenched approximation in lattice QCD 
study [42]. Especially the classical supergravity descrip- 
tion of D4 branes is tractable, followed by the local ap- 
proximation of the strings and also the suppressions of 
the string loop effect. These conditions in the gravita- 
tional side around the D branes give the constraints for 
the gauge theory side as the QCD on the surface of D 
branes as 



which is achieved by qym —>■ 0, Nc 00, and 'tHooft 
coupling: A = gy^Nc fixed and large. In this sense, 
the strong-coupling large-A^c QCD is found to be dual 
with the classical supergravity of D4 branes with probe 
D8 (D8) branes, which is one of the realization of the 
strong-weak duality between the gauge theory and grav- 
itational theory through the holography. Therefore, by 
analyzing the effective action of D8 branes with D4 su- 
pergravity background around the D branes, one can an- 
alyze the nonperturbative aspects of QCD on the surface 
of D branes. 

Then we start formal discussions from the Nf = 2 
non-Abelian Dirac-Born-Infeld (DBI) action of probe D8 
brane with D4 supergravity background as a probe ap- 
proximation. After dimensional reductions, the nine- 
dimensional DBI action of the probe D8 brane with 
the D4 supergravity background can be reduced into a 
five-dimensional Yang-Mills theory, belonging to the flat 
four-dimensional Euclidean space-time x, i.e., ^0,^3 and 
the other fifth dimension z with curved measures as fol- 
lows [3]: 

+K{z)F^,F^,} + 0{F^), (2) 

where Aj^i is the gauge field and Fmn — QmAn — 
On Am + i[AM,AN] {M,N = - 3,z) is the field 
strength tensor in five-dimensional space-time {xo,^3,z) 
of the probe D8 branes. In the action (2), Mkk = 1 unit 
is taken, and the overall factor k is defined as 



(Note here that we use the value of k as a half of that in 
Rcf. [9], taking away the misleading factor 2 in Eq.(16) 
of Ref . [9] . All the formula and numerical results can be 
scaled by the factor k, so that the discussions in Ref. [9] 
are not altered.) The functional K{z) = 1 -t- in the 
action (2) expresses the nontrivial curvature in the extra 
fifth dimension z induced by the supergravity background 
of the D4 brane. The gravitational energy of D8 brane, 
i.e., SEs^Iam-^o is subtracted in the action (2) as the 
vacuum relative to the gauge sectors. 

In the D4/D8/D8 multi-D brane configurations, color 
quantum number is carried only by the iVc-folded D4 
branes. Therefore, after the supergravity description of 
D4 branes, there is no colored particles from the fluc- 
tuation modes of open strings on the residual probe D8 
branes. This is regarded as some holographic manifesta- 
tion of "color confinement" in low-energy scale of QCD. 
In fact, gauge field Am=o~3.z in the action (2) is color- 
singlet and obeys the adjoint representation of U{Nf) 
group, eventually producing the meson degrees of free- 
dom after some proper mode expansions in the holo- 
graphic QCD. 

In this paper, we treat the non-trivial leading order 
of and 1/A expansions in the holographic QCD as 
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the five-dimensional Yang-Mills action (2) with 0{F'^). 
In general, there also exists the Chern-Simons (CS) term 
to avoid anomalies in the superstring theory. By intro- 
ducing u! meson degrees of freedom as the U{1) sector 
of p meson fields in the holographic approach, lu meson 
indirectly couples with pions via (axial) vector mesons 
after the proper mode expansions of CS term, which is 
regarded as a general representation of the GcU-Mann - 
Sharp - Wagner (GSW) model [7]. The roles of cu mesons 
for low-energy meson dynamics and also chiral solitons 
have been traditionally examined in some of QCD phe- 
nomenologies [43]. Actually, however, the CS term is 
found to be 0{X'^), i.e., the higher order contributions of 
'tHooft coupling expansion relative to the 0{F'^) of the 
Yang-Mills action (2) with O(A^), which is manifestly 
shown in holographic QCD. Furthermore, most part of 
CS term includes one time-derivative of pion fields and 
docs not affect the static properties of hedgehog solitons. 
Therefore we neglect the CS term along the discussions 
below for the argument of the non-perturbative (strong 
couphng) properties of QCD. 

In the holographic approach, the pion field is intro- 
duced as the Wilson line of the fifth gauge field A^, i.e., 
a path-ordered product of the fifth gauge field along the 
z direction [3, 9, 31] as 

U{x'') ^ Pcxpl^-i J dz'A^(a;^,z')| € f7(iV/). (4) 

One can also introduce the variables ^±(x^) as 



can be mode-expanded into the four-dimensional parity 
and G-parity eigen states with proper complete orthogo- 
nal basis 4'±{z) and tpniz) (n = 1, 2, • • • ) as follows [3, 9]: 

+ Y,B(-\x.)Mz), (7) 



1 



(8) 
(9) 



where and are left and right currents of pion fields, 
respectively. The basis '4'±iz) are introduced to support 
whole of the gauge field Afj_{xN) at the boundary z 

±oo as 'ip±{z ±oo) ~ 1 and 'ijj±{z Too) = as 



"00 (z) = — arctanz. 



(10) 
(11) 



In order to diagonalize the five-dimensional Yang-Mills 
action (2)with the induced measures K{z)~^^^ and K{z) 
in the fifth dimension z, the basis ipn ('^ = Ij 2, • • • ) are 
taken to be the normalizable cigcn-function satisfying 



-A-(z)V3^ Ik{z)^ } = A„Vn, (Ai < A2 < • • • ) 



(12) 



e±i(a-^) = Pcxp|-j 



dz'A,{x^,z') \ e U{Nf), 



(5) 



with normalization condition as 

K I dzK{z)~'^/^1prntpn = 5„ 



(13) 



where zo(x^j) is a single- valued arbitrary function of a;^ 
Then the pion field (4) can be written as 



U{x^^)^£,:,\x^)^.{x^), 



(6) 



some resemble formula of which can also be found in the 
traditional approach of hidden local symmetry [4] . 

Now we take "A^ = gauge" and also "^^^(a;^) 
^_(a;^)(= i{x^)) gauge" for the U{Nf) gauge symmetry 
in the action (2) of the probe D8 brane. The A^ = 
gauge is similar to the unitary gauge in the non-Abelian 
Higgs theory; fifth gauge field A^ performs as a scalar 
field in four-dimensional space-time Xp=o~3 7 and it is 
eaten by the four dimensional gauge field to give the 
mass generation of gauge field, especially the (axial) vec- 
tor mesons as a part of A^^. In this sense the masses of 
the (axial) vector mesons come from the Higgs mech- 
anism in five-dimensional space-time with the U(Nf) 
gauge symmetry breaking on the probe D8 brane. The 
i+^ix/j,) = £,-{Xfj,) gauge is also essential to get the low- 
energy effective theory of QCD with proper parity and 
G-parity classification in a manifest way [9] . With these 
gauge fixings, the five-dimensional gauge field Af^{x]y) 



In the holographic model, the fields ■* in 

the mode expansion (7) arc regarded as (axial) vector 
mesons, belonging to the adjoint representation of the 
U{Nf) gauge group as = B^T"" . By substituting 

the expansion (7) into the action (2), the mass of B^"-* 
field is found with the eigenvalue of oscillating fifth ba- 



sis in Eq.(12) as 



A„, indicating that the origin 



of meson mass is the oscillation of meson wave func- 
tion in the extra fifth dimension. Furthermore, A^{xm) 
is the five-dimensional vector and ipn{z) are the parity 
eigen state in the z direction as ^n{—z) = {—)"'~^'ipn{z)- 
Therefore, from the mode expansion (7), Bjr^ fields have 
four-dimensional parity transformation as B^\^Xi,) — > 

{—y^Blr\xi,). These consideration indicates that vec- 
tor and axial vector mesons appear alternately in the 
excitation spectra about index n as B^}^ = p^^B^if' = 



n' B^'^' = a' B^^) = n" 



In our study, we construct the four-dimensional meson 
effective action only with pion field U{xi,) and p meson 

field b'iP {x^) = p^{x^) below the Kaluza-Klcin mass 
scale Mkk ~ IGeV. Recall that there appear infinite 
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number of non-QCD modes with large mass 0{Mkk) 
in holographic QCD, e.g., scalar fields, gluino, and also 
the Kaluza-Klein modes, discussed in the first part of this 
section. Therefore, the D4/D8/D8 multi-D brane system 
can be viewed as QCD as far as low energy phenomenol- 
ogy below ~ Mkk is concerned. In this sense, Mkk 
plays the roles as the ultraviolet cutoff scale of the the- 
ory, so that we include the meson degrees of freedom be- 
low Mkk for the baryon analysis in later sections, called 
"truncated-resonance model" . 

Actually, the appearance of Mkk scale with finite value 
seems to be essential in the recent holographic analysis. 
In the framework of AdS/CFT correspondence without 
Kaluza-Klein compactification, the = 4 SUSY and its 
resulting conformal symmetry protect the emergence of 
the dimensional quantities like string tension and chiral 
condensate at the ground state, because of the cancel- 
lation of the radiative corrections between bosons and 
fermions. In this sense, SUSY breaking is at lest needed 
to be dual of QCD with confinement and chiral symme- 
try breaking as its vacuum nature. In the Sakai-Sugimoto 
model, Nc sheets of D4 branes are Kaluza-Klein compact- 
ified with radius Mj^^, and the SUSY is completely bro- 
ken by the field boundary condition along the compacti- 
fied direction. Therefore confinement and chiral symme- 
try breaking could occur as the non-SUSY gauge theory. 
In fact, the compactificd D4 brane is "non-BPS", hav- 
ing a "horizon" in the supergravity description. In the 
holographic framework, confinement and chiral symme- 
try breaking do occur on this horizon, giving the lost of 
"colored" information, and also the geometrical connec- 
tion of D8 and D8 branes. As a whole, the appearance 
of Mkk would be essential for the holographic model to 



be dual of realistic QCD. 

Now, one may regret about the finiteness of Mkk as 
almost IGeV, which is comparable with the QCD mass 
scale Aqcd • In the holographic approach, the QCD mass 
scale is introduced by the experimental inputs for the 
pion decay constant and p meson mass as = 92.4MeV 
and nip = 776.0MeV. These experimental values come 
from our hadronic world with Nc = 3 and large but fi- 
nite 'tHooft coupling A, which may infringe the condition 
(1) to give the effects of string length and string loops in 
the gravitational side. Therefore, too large Mkk can not 
be taken to neglect the internal structure of strings on 
the surface of D4 branes compactificd with radius M^^, 
which might eventually gives the scale A/kk ^ IGcV. 
These considerations indicate that, by including the ef- 
fects of string length and loops in the gravitational side, 
Mkk would be taken sufficiently large relative to Aqcd 
with fixed and nip. Anyway, Mkk ~ IGeV essen- 
tially appears as the ultraviolet cutoff scale in the sys- 
tem of probe D8 brane with D4 "classical" supergravity 
background to be dual of QCD with proper dimensional 
quantities. 

By neglecting the higher mass excitation modes of (ax- 
ial) vector mesons rather than p meson sector with n = 1 
in the expansions (7), the five-dimensional gauge field 
Afj,{xN) can be written as 

A^^Ixn) = lp{x„)tjj+{z) + rp{xv)ip-{z) -)- pp{xv)ipi{z). 

(14) 

By taking this mode expansion (14) with the — Q 
gauge, five-dimensional field strength and F^p can 
be written as 



Fp,, = dpA^ - dyAp -t- i \Ap, A^\ 

= {dplu - djp) V'+ + {dpri, ~ di,rp) -0- -f {dpPu - d^Pp) tpi +i{[lp, + [rp, rv]ip'^ 

+i{{[lp,ri,] + [rp,l^])ip+tl)^ + {[lp,p,.] + [pp,l^]) i'+tpi + {[rp,p„] + [pp, r„]) tp-ipi} 

= -i[ap, av]ip+ip- + [dppu - d^Pp) ipi + i[pp, Pujtpl + i{i[atJ., Pu] + [pp., a^]) ipo^'i 
+ i[Pp,p,] + [pp,(3,])^i}, 
Fzp = dzAp = UpdztpQ + PpdzTpi, 



(15) 
(16) 



with axial vector current ap and vector current Pp of 
pion field as 



ap{x^) 

Pt,{Xy) 



'rrp{x^)}. 



(17) 
(18) 



with 0{F^), we eventually get the four-dimensional Eu- 
clidean meson effective action with pions and p mesons 
from holographic QCD as follows (derivations in more 
detail can be found in our previous paper [9]): 

c _ qDBI cDBIi 

•Jeff = ^DS ^ '-'D8 \Am-*0 



In the derivation of (15) and (16), we have used the = k f d'^xdztTl-K{z)^^^^Fpi,Fpi, + K{z)FpzFpz\ 

Maurer-Cartan equations, — -I- Z^] = and J i J 

dnr,, — duVn + i\ru,r^] = 0. By substituting Eqs.(15) (19) 



dpTi, - d^Vp + i[rp,r^] = 0. By substituting Eqs.(15) 
and (16) into the fivc-dimensional Yang-Mills action (2) 
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J TT 

4 „ 

+32 
+«.94 

1 

-77.96 



32e2 



^ d''.Ttr{(5pP,. - - -g4p J d'^xti- [p^,pvf - igi j d'^xtr {[a^.ay] {d^j,pu - duP^^)} 

d'*i-tr{[a^,Q;^] +.93 j d'^xtr{[a^,au\ + [p^n 

d'^xiY{[a^j,,pi,] + [p^,,a^]f -]^g^ j d'^xiv {[13^, p^,] + [p^,, I3u\f , 



(20) 



where is the 1-form of pion fields as 

L^{x,) = -U\x,)d^U{x,). 



(21) 



There exist twelve kinds of coupling constants: /tt, mp, e, 
53p? .94p s-nd 5i^7 in the action (20). However all the cou- 
pling constants are uniquely determined from the proper- 
tics of meson wave functions in the extra fifth dimension 
z, i.e.; the complete orthogonal basis i!±{z) and i!i{z) 
with oscillating eigen value Ai as follows: 



J TT 

4 



,2 _ 



1 



16e2 

53p 

91 

92 
93 
94 
95 
96 

97 



K J dzK{z 
m\ = Ai, 
K j dzK{z 

dzK{z 



(dA^f = -, (22) 

TT 

(23) 
(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

52, (32) 

(33) 

The holographic model has two parameters k(= ^j^fj) 
and the Kaluza-Klcin mass Mkk as the ultraviolet cutoff 
scale of this theory, k appears in front of the effective ac- 
tion (2) because the effective action of D8 branc with D4 



dzK{z 
K J dzK{z 
K J dzK{z 
K J dzK{z 
K J dzK{z 
K J dzK{z 
K J dzK{z 
K / dzK{z 



-V3^^(l-^+) 

-'^H'1=93P, 
= 1. 



Si, 



supergravity background expanded up to 0{F'^) corre- 
sponds to the leading order of 1/Nc and 1/A expansions. 
Therefore, by fixing two parameters k and Mkk to adjust 
experimental inputs for f^^ and rup, then all the coupling 
constants (22)^(33) are uniquely determined through the 
background of extra fifth dimension z. Such uniqueness 
of the action is one of the remarkable consequences in 
holographic QCD. 



III. BRANE-INDUCED SKYRMION ON 
AND 

In this section, we discuss baryons and baryonic matter 
in holographic QCD. In Sec. Ill A, we describe the baryon 
as a chiral soliton in the four-dimensional meson effec- 
tive action in (20) including pion and p meson fields 
derived from holographic QCD. We call this topological 
soliton as the "brane-induced Skyrmion" . The hedgehog 
mass of the brane-induced Skyrmion on the flat coordi- 
nate space R'^ is derived, which is originally given in our 
previous paper [9]. In Sec. IIIB, we newly discuss the 
baryonic matter in holographic QCD by analyzing the 
system of single brane-induced Skyrmion on the three- 
dimensional closed manifold S'^. Through the projection 
procedure from the flat space R'^ onto the curved space 
the hedgehog mass and the Euler-Lagrange equations 
of the brane-induced Skyrmion on are derived. All 
the numerical results and their physical interpretations 
are presented in Sees. IV, V, VI. 



A. Brane-induced Skyrmion on R'^ 

In this work, we describe the baryon as the four- 
dimensional chiral soliton, i.e., the Skyrmion in holo- 
graphic QCD. To sec the validity of this approach for 
the baryon, we now compare the meson effective action 
induced by holographic QCD with that in a chiral pertur- 
bation theory (ChPT) as a low-energy effective theory of 
QCD [44]. The ChPT is phenomenologically constructed 
respecting the symmetry of QCD, the chiral symmetry 
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and the Lorentz invariance in the four-dimensional space- 
time. With these symmetry constraints, there are three 
possible terms as the four-derivative terms of pion fields: 

tr[L^,L,f, tr{L^,L,f, trid^L.f, (34) 

where = jU^dfj^U is the 1-form pion fields in Eq. (21). 
The first term tr[L^, L^J^, called the "Skyrme term" [45], 
is to give the stability of the Skyrme soliton solution with 
finite size in the coordinate space. On the other hand, 
the other two terms are known to give the instability 
of Skyrme solitons [46]. The symmetry constraints in 
the CliPT can not determine which terms should appear 
because all the terms in (34) are chiral symmetric and 
Lorentz invariant. Therefore Skyrme deliberately takes 
only the first term tr [L^,Ly]^ in the meson effective ac- 
tion as an effective "model" for the baryon as the chiral 
soliton, which was called the "Skyrme model" [45]. 

Now, by starting from the holographic QCD with the 
five-dimensional Yang-Mills action (2) of the probe D8 
brane, one can find only the Skyrme term without the 
other two in (34) , which is manifestly seen in the action 
(20). Actually the five-dimensional Yang-Mills action (2) 
with 0{F^) includes two-time derivatives at most, so that 
the appearances of the other two terms: tr{L^, L^}^ and 
tr{df^L^)^ with four-time derivatives are forbidden at the 
leading order of 1/Nc and 1/A expansion in holographic 
QCD. These comparisons with the chiral perturbation 
theory clearly indicate that holographic QCD is not just 
the low-energy effective theory of QCD only with the 
constraint of symmetries in four-dimensional space-time: 
actually, it obeys the U{Nf) symmetry extending to the 



extra fifth dimension z. Furthermore, one can see that 
the chiral soliton picture for the baryon is now supported 
by the holographic approach, retaining the direct connec- 
tion with QCD. With these considerations, we employ the 
concept of chiral soliton picture for the baryon analysis 
in holographic QCD [9]. 

Now we begin with the hedgehog Ansatz for pion field 
U (x) and p meson field (x) as a baryon configura- 
tion [9]: 

U^i^) = e-»^"^W, (xa = ^,r= |x|) (35) 

{G{r)^G{r)lr) (36) 

where is Pauli matrix, and F{r) is a dimensionless 
profile function of the pion field with boundary conditions 
F(0) = TT and F{oo) = 0, giving topological charge equal 
to unity. Ansatz (35) means 7ra(x) = XaF{r) for the pion 
field. G(r) is also a dimensionless profile function of the 
p meson field. This Ansatz for p meson field is also called 
"Wu-Yang-'tHooft-Polyakov Ansatz" [47], and the same 
configuration Ansatz can be seen for the gauge field of 
the 'tHooft-Polyakov monopole [48] . 

By substituting the configuration Ansatz (35) and (36) 
into the four-dimensional meson effective action S'off in 
(20) with the Euclidean metric, we can get the static 
hedgehog mass of a brane-induced Skyrmion on the flat 
space R'^ as follows (detailed derivations can be found in 
our previous paper [9]): 



E[F{r),G{r)] 



rc-DBI _ cDBI 



D8 I^A/-»Oj[^j,jgpi^Qg 



ATTdrr^ ■ e[F{r),G{r)], 



•I TT 



. £[F(r), G(r)] = ^ [2 (r^F'^ + 2 sin^ F)] 



2 

53 



\[?>G^ + 2rG{G')+r^G'^^ 



93p 



sinF- (G + rG'] + sin^ F • 



f6 



Wsm^F ■ (1 -cosF) G/r 



32e2 
f6rG^ 

G/r} 



16sm^ F I 2F'2 
f 



sin^F 



T54p 



f6r^G" 



52 



16sin^F- G" 



.94 



f6(f -cosF) G^ 



55 



f6r(f -cosF) G^ 



56 Wr^F'^G'^ +57 8(1-cosF)^G2 , 



(37) 



(38) 



where F' ^ ^) and G' ^ ^ 

dr ^ or ' ar 



dGjr) 



Now we take the "Adkins-Nappi-Witten (ANW) unit" for energy and length as Fanw = It ^^'^ ^anw = 7^ [25], 
and we rewrite all variables in this ANW unit as F = -r^ — F and r = — - — r. By taking this scaled unit, the 

-C/ANW rANW .JO I 

hedgehog energy density (38) of single brane-induced Skyrmion on R'^ can be rewritten as follows (overlines of F and 
r below are abbreviated for simplicity) : 



£[F(r),G(r 



(r^F'^ + 2sin^F) + 2 ( -A 



4r^G^ 



sin^ F f 2F'2 



sin^ F 
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(26^) i [S {SG^ ^ 2rG(G') + r^G'^}] - {2e^) g^p 



16rG' 



(2e') ^ff4p 



ler^G^ 



(2e2) .91 



16 {f' sin • (g + rG') + sin^ F ■ G/r } - (2e^) 52 
IGsin^i^- (1 - cos F) G/r] - (26^) .94 [l6(l-cosF)G 



16 sin^ F ■ G' 

■2 



2e^) 55 16r (1 - cos F) G^ + (2e^) 56 



ler^F'^G^ 



2e^ 57 



8(1-cosF)^G2 



(39) 



Here we comment about a scaling property of the 
brane-induced Skyrmion. The holographic QCD has just 



two parameters: k(= 2iq^) and Mkk, so that the pion 
decay constant /jr, the p meson mass nip, and the Skyrme 
parameter e in Eqs. (22)^(24) can be explicitly written by 
K and Mkk in holographic QCD as 



2^1 -Mkk, 



-1/2 



(40) 

(41) 
1 



4^/(^57 V^' 
(42) 



where the energy unit Mkk is recovered. By using these 
relations (40)'--'(42), ANW unit for energy and length, 
i.e., Fanw and tanw can be written by k and Mkk as 



U 

Fanw = 7;- = const • kMkk, 
2e 

1 1 
'"ANW = —r = const • — — . 
eU Mkk 



(43) 
(44) 



The five-dimensional Yang-Mills action (2) with 0{F^) is 
proportional to k(= 2i^:i ) , as the leading order of 1/Nc 
and 1/A expansion. Furthermore Mkk is the sole en- 
ergy scale of the holographic approach. Therefore, in the 
energy unit Fanw(oc kMkk), the total energy appears 
as a scale invariant variable. In fact, by introducing the 
rescaled p meson field G(r) as 



G(r) ^ ^G(r), 



(45) 



and considering the K-dependence of the basis ipi as 

-01 oc -7= in the normalization condition (13), one can 
v ^ 

analytically show that every energy density in each term 
of Eq.(39) and meson field configurations F{r) and G(r) 
are scale invariant variables, being independent of the 
holographic two parameters, k and Mkk- 

With the considerations above, we give most discus- 
sions below in the ANW unit as the universal features 
of baryonic matter in holographic QCD, being indepen- 
dent of the definite values of and nip. The recovering 
of physical unit with the experimental inputs for and 
nip is discussed in Sec. VI, with respect to the critical 
densities of the phase transitions in the baryonic matter 
within holographic approach. 



B. Brane-induced Skyrmion on S 

Now we study the baryonic matter in holographic QCD 
by analyzing the system of single brane-induce Skyrmion 
on a three-dimensional closed manifold S^. 

In this study, we consider the baryonic matter with 
large Nc because holographic QCD is the large Nc effec- 
tive theory, derived from the classical supergravity justi- 
fied in the large Nc and large 'tHooft couphng [49]. Ac- 
cording to the general analysis of large- iVc QCD, a baryon 
mass is found to become 0{Nc) [8, 50], so that its kinetic 
energy becomes 0{N~^). As for the quantum effects of 
the baryonic matter, zero point quantum fluctuation en- 
ergy Eq and baryon mass splitting Am in the isospin 
projection correspond to the higher-order effects of 1/A^c 
expansion: Eq ^ 0{N°) and Am - 0{N-^) [25]. There- 
fore, with large-TVc condition, we can consider that the 
kinetic energy and quantum effects are suppressed rela- 
tive to the static mass, and the baryonic matter comes 
into the "static Skyrme matter" . 

Such static Skyrme matter was first analyzed by Kle- 
banov [26], placing Skyrme soliton configurations peri- 
odically along the three-dimensional cubic lattice, which 
could be related with "nuclear crystal" with pion conden- 
sation in the deep interior of neutron stars. Therefore, 
by analyzing the static Skyrme matter, one can see some 
typical features of baryonic matter with large- A^c condi- 
tions. 

In this paper, we take certain mathematical trick to 
analyze such static Skyrme matter suggested by Man- 
ton and Ruback [27]. To analyze the features of multi- 
Skyrmion system on the flat coordinate space R'^. they 
alternately treat the system of single Skyrmion on a 
three-dimensional closed manifold with finite radius 
R as shown in Fig. 1. Actually, the multi-Skyrmion sys- 
tem on R'^ and the system of single Skyrmion on 
can be related with each other through the compactifi- 
cation of the boundary for a unit cell on R^ shared by 
the single Skyrmion as in Fig. 1. The interaction be- 
tween the baryons in the medium on R"^ is simulated by 
the curvature of the closed manifold S^. The baryon- 
number density can be represented as = l/27r^i?'^ 
on S^, so that, as the size of decreases, the increase 
of the baryon-number density in the medium is repre- 
sented in this modeling. Actually, as the radius R of 
decreases, the energy density of single baryon is found 
to delocalize due to the medium effects in the baryonic 
matter [27], and, below the critical radius i?crit of S^, the 
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Skyrme matter 




R <Rcnt 

Pb > Pb(c) 



R '. small 

i 

Pb : large 



compactification 




Identity map 



FIG. 1: Schematic figure of the static Skyrme matter on a flat 
coordinate space R^, and the system of single Skyrmion on a 
closed manifold with finite radius R. The static Skyrme 
matter on Yi? and the system of single Skyrmion on can be 
related with each other through the compactification of the 
boundary for a unit cell on Y{? shared by the single Skyrmion. 
The decrease of the radius R of S"^ represents the increase of 
the baryon- number density Pb{= {2tv^ R^)~^) in the medium 
in this modeling. For R < Rait as a critical radius, i.e., 
Pb > Pb(c){= {27r^(7icrit)^}~^) as a critical density, the en- 
ergy density of the single Skyrmion becomes uniform distri- 
bution as the "identity map" discussed in Eq. (67), which is 
called the "delocalization phase transition" . 



energy density of the baryon coincides with the uniform 
distribution as the "identity map" , which is called the 
"delocalization phase transition" shown in Fig. 1. Such 
delocalization phase transition in the Skyrme model can 
be related with the deconfinement of the baryon and also 
the chiral symmetry restoration in QCD, which will be 
inclusively discussed by taking the order parameters in 
Sec. V. With these considerations, by analyzing the sys- 
tem of single brane- induced Skyrmion on S^, we can see 
some typical features of baryonic matter in holographic 
QCD. Especially, by comparing the standard Skyrmion 
without p mesons and the brane- induced Skyrmion on S^, 
the roles of (axial) vector mesons in high density phase 
of baryonic matter can be discussed from the holographic 
point of view. 

Now we introduce the projection procedure from the 
fiat space onto the curved space 5*^ [51], to get 
a hedgehog mass of a brane- induced Skyrmion on S^. 
First, we consider the three-dimensional orthogonal space 



R in polar coordinates as 



X = {z,x,y) 

= {r cos 6, r sm9 cos (j),r sin9 sine 

= (r, 0, 0)3dim. polar- 



(46) 



The integral operator dx and the derivation d can be 
written in polar coordinates as 



dx = (dr,d6',d0)3dim. polar 

= {dr, rde, r sin 6'd0)3dim.poia 
dr r do 



(47) 



d 



r sin dcj) 

so that the differential operator d can be written as 

/ 3dim. polar 



(48) 



d 



d_ ld_ 

dr^ r dO^ r sin ( 



(49) 



Second, we consider the four-dimensional orthogonal 
space R* in polar coordinates as 

X = {t,z,x,y) 

= (i?cos 8, i?sin O cos 6*, sin 8 sin cos 0, 
R sin 8 sin 6 sin 0) 

= (i?, 8, 6*, 0)4dim.polar- (50) 

The integral operator dx and the derivation d can be 
written in polar coordinates as 



dx = 



{dR, dQ, dO, (i(/))4dim. polar 

{dR, RdQ, R sin 8^6*, R sin 8 sin 



4dim. polar ; 

(51) 



,-9 ,^1 d 



dO 



1 



d 



i? sin 8 961 



1 



d 



R sin 8 sin 6 dcj) ' 
so that the differential operator d can be written as 



(52) 



d 



1 



d 



1 



d 1 d 



dR Rde RsinOde i?sin8sin6i, 



d 



4dim. polar 



(53) 



Now, by limiting the four-dimensional orthogonal 
space R'* onto the surface of a three-dimensional closed 
manifold with fixed radius R, the coordinate t in 
Eq. (50) becomes dependent on the other coordinates 
{z,x,y). Furthermore, dR and ^ can be regarded as 
zero in Eqs. (51) and (53) because the radial coordinate 
R is fixed on S^. Therefore, by comparing Eqs. (46), 
(47), (49) on R^ , and Eqs. (50), (51), (53) on with 
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fixed radius R, wc find the projection procedure from 
to S'^ as follows: 



r 
dr 
d_ 
dr 



i?sin6, 
Rde, 

Rde' 



(54) 
(55) 

(56) 



Recall that the hedgehog mass on with energy den- 
sity in Eq. (39) can be written with its explicit arguments 



for the energy density as 



'o 



F(r),G(r),|-F(r),|-G(r),r 



(57) 



By applying the projection procedure (54), (55) and (56) 
to the hedgehog mass on R"^ in Eq. (57), we can get the 
hedgehog mass on as 



E 



ATTRdQR^ sin^ 9 • e 

AirdrR^ sin^ — • e 
R 



F(R sin 9), G(R sinQ), -^F{R sin 6), -^G(R sin 6), R sin 9 
R o(3 R a9 

Fir),G{r), — F(r), — G(r), R sin ^ 
or or R 



In Eq. (59), we introduce a new variable r as the arc length on as 

r = i?9, 



(58) 
(59) 

(60) 



and r-dependent dimensionless functions i^(i? sin -^) and G(i? sin -^) are renamed again as F{r) and G{r). Therefore, 
by comparing Eq. (57) on and Eq. (59) on S'^ , we can get the simple projection procedure for the hedgehog energy 
density from R'^ to S*'^ as 



drr 



F{r),G{r),-^F(r),-^G{r),r 
or or 



R 



F(r),G(r),— F(r),—G(r),i? sin 



where the topological boundary for the chiral field F{r) in (35) is also projected on 5*^ as 

F{0) = TT, F{ttR) = 0. 



(61) 



(62) 



Now, by applying the projection procedure (61) for the hedgehog energy density on R^ in Eq. (39), we can eventually 
get the hedgehog energy density on with ANW units as follows (Note here that the dimensional profile function 
G(r) = G{r)/r on R'^ is to be naturally introduced on through the projection procedure as G(r) = J'-^K ) : 



£;[F(r),G(r)] 



AndrR sin 



R 



r[F(r),G(r)], 



RHin'L.e[F{r),G{r)] = (r^ sin' ^ ■ F'^ + 2sm^ F] + 2 [ ^ 

R \ R ^ \ Jn 



4i?2 sin^ - • G^ 
R 



sin^ F 2F'- 



sin^F 
i?2 sin^ ^ 



2e^)^ 



8 I f 2 + cos^ ^) G^ + 2i? sin ^ cos ^ • G(G') + i?^ sin^ ^ ■ G'^] 
LV R/ R R R i 



2e') 53p 



R 

16i?sin- • G^ 
R 



16R' sin^ - • G^ 



2e^) 51 lejF'sinF- (^cos • G + i? sin • G') + sin^ F • G/ (^i? sin -^^ | 

16sin2F • (1 - cosF)G/ (^Rsin^^ 



2e') 92 



IGsin^F- G' 



R 



2e^).g4 16(1-cosF)G2 +(26^)^5 16i?sin - • (1 - cosF) G^ 
L J L it 



2e2 56 



16E^ sin^ 



R 



26^) .97 



8(1-cosF)^G2 



(63) 



(64) 
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where F' = ^^^7^(= ^^7^) and G' = ^^^7^(= ^^7^)- Here measures i?^ sin^ newly appear in the energy density 
(64) relative to Eq.(39) on the flat space R^, indicating the existence of the curvature of the closed manifold 5''^. We 
also construct the Euler-Lagrange equations for the pion field F{r) and p meson field G{r) from the energy density 
in Eq. (64) as follows: 



5E 



d f 6E 



An {SF{r) dr \SF'{r) 

-4i?sin — cos — • F' - 2i?2 sin^ — • F" + 4sini^ • cosi^ 
R R R 

I -4 sin F • cos F • - 4 sin^ F • F" + 4 sin^ F • cos F/ [r^ sin^ ^) | 



(2e2) gi 

(2e') 92 
(2e2) .94 
(2e2) .96 
(2e2; 



16 <^ 2sinF • cosF • G/ LRsin 



) 92 [16 (: 



16 2sinF- cos F-G^ 



R 



- (2e') 53 



sinF ■ 



- sin - • G + 2cos - • G" + i? sin - • G" 
R R R R 



16 (sinF + 2sinF • cosF - 3sinF • cos^ F) G/ (i?sin — ) 



16 sin F • G' 



26^)55 16(^i?si 



sin — sin F • G^ 
R 



16 



-4i?sin - cos - • F'& - 2i?2 sin^ - • F"& - AR^ sin^ - • F'GG" 
R R R R 



).97 8 |2(l-cosF)sinF-G2| 



0, 



RJ 



(65) 
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Now one can show that the Euler-Lagrange equations 
(65) and (66) always have the analytical solution for ar- 
bitrary value of radius R of S'^ as 



F(r) =7r--, G(r) 



0, 



(67) 



which is the "identity map" solution in the standard 
Skyrme model without p meson field [27] . In the Skyrme 
model, the energy density of the baryon is generated by 
the spatial gradient of the pion fields. Therefore, such 
linear configuration as in Eq.(67) gives the uniform en- 
ergy distributions with its classical mass as 



which can be given by substituting the linear config- 
uration (67) into the energy density (64). Hence the 
change of absolute minimum solution from the localized 
Skyrmion into the identity map (67) can be regarded as 
a signal of transitions from localized phase to the uni- 
form phase of the baryonic matter, referred as the delo- 
calization phase transition. Note that the identity map 
solution (67) has no p meson configuration: G(r) = 0. 
Therefore, if the identity map is realized as the absolute 
minimum solution, it indicates that p meson field abso- 
lutely disappear in high density phase, which is discussed 
with numerical results in Sec. IV. 



Fid = ( i? + I 6^^ 



(68) 



13 



IV. BARYON NATURE IN DENSE QCD 

In this section, we show the numerical resuhs and 
discussions about the baryon nature in dense QCD by 
solving the Euler-Lagrange equations (65) and (66) de- 
rived from holographic QCD. The energy density and the 
field configuration profiles of the baryon arc analyzed in 
Sec. IV A. A new striking picture of "pion dominance" 
near the critical density is proposed in Sec. IV B. The 
mass and root-mean-square mass radius of the baryon 
are examined in Sec. IV C. We explain the "swelling" 
mechanism of the baryon in the general context of QCD 
in Sec. IV D. Through all of the sections below, by com- 
paring the Skyrmc model without p meson field and the 
brane-induced Skyrme (BIS) model on 5*^, the roles of 
(axial) vector mesons in the baryonic matter are dis- 
cussed from a holographic point of view. 

A. Energy density and field configuration profiles 

In this section, we discuss the baryon-numbcr density 
dependence of the energy density and the field config- 
uration profiles of single baryon for the Skyrme model 
and the BIS model, by changing the size of the manifold 
S'^ . Actually, the baryon- number density is represented 
as pb = on 5*^, so that, as the radius R of 

decreases, the increase of total baryon-numbcr density in 
the medium is represented in this modeling. 

First, we show in Fig. 2 the energy density of the baryon 
for the Skyrme model with radius R oi . One can see 
that the energy density tends to delocalize as R decreases, 
which can be regarded as some medium effects in the 
dense baryonic matter. Below the critical radius, R < 
-^crit'^™'' = a/2, the energy density distribution of the 
baryon exactly coincides with the uniform one as identity 
map denoted by the dashed lines in Fig. 2. The transition 
from localized energy density into uniform one is called 
the "dclocalization phase transition" . 

Next we show in Fig. 3 the same plot in the case of the 
BIS model. The energy density of the baryon also tends 
to delocalize as R decreases. However the delocalization 
along with the decrease of R is delayed relative to the 
Skyrme model without p meson field, and the energy dis- 
tribution is still localized even around R - R^^^^^ = V2, 
which is the critical radius for the Skyrme model. In 
fact, the heavy p meson field appearing in the core re- 
gion of the baryon is to provide the attraction with the 
pion field, which leads to the shrinkage of the total size of 
the baryon [9]. Therefore, the smaller radius of S^, i.e., 
the larger baryon-number density is needed for the BIS 
model to give the delocalization phase transition, which 
is discussed in Sec. VI with recovering the physical units. 
Then, below the critical radius, R < Rf^^ = 1.19, the en- 
ergy density distribution of the baryon in the BIS model 
becomes the uniform one as identity map (67). These 
delocalization phase transitions in the Skyrmion picture 
can be related with the dcconfincmcnt of the baryon and 



also the chiral symmetry restoration, which will be dis- 
cussed with the order parameters in Sec. V. 

We can also show in Fig. 4 the comparison between the 
energy density of the baryon and p meson contributions 
in the interaction terms of the BIS model (64) for R ~ 
4.0, a/2, and 1.19. For R = 4.0 and a/2, one can see the 
manifest contributions from the p meson field in the core 
region of the baryon. On the other hand, at the critical 
radius R = R^^^ = 1.19, all the contributions from the p 
meson field absolutely disappear in the uniform phase. 

We also show the field configuration profiles of the pion 
and the p meson in Fig. 5 and Fig. 6, respectively, in the 
BIS model for R = 4.0, a/2, and 1.19. As R decreases, 
the pion field F{r) in Fig. 5 approaches to the linear 
configuration as the identity map in Eq.(67), giving the 
uniform energy density distribution. As for the p meson 
field in Fig 6, the amplitude of p meson field G(r) tends 
to decrease as R decreases, and it absolutely disappears 
below the critical radius, R < Rf^^ = 1.19 as the identity 
map in Eq. (67). These results indicate that the ampli- 
tude of the p meson field decreases as the baryon-number 
density increases in the medium, and it absolutely disap- 
pears and only the pion field survives near the critical 
density. 

B. Pion dominance near critical density 

In the previous section, we find that the p meson field 
would disappear near the critical density. Now we pro- 
pose a conjecture that such disappearance of the p me- 
son field near the critical density can be generalized to 
all the other (axial) vector meson fields: ai, p' , a'l, p" ■ ■ ■ , 

denoted by the field B^ff'{xv) in Eq. (7) by the following 
reasons within holographic approach: 

1) The kinetic term of (axial) vector meson field B^"'' 
on the closed manifold is proportional to R~^ 
as 

tr {a^B^") - cx (69) 

which can be derived from a simple dimensional 
analysis. In general, the kinetic energy indicates 
the "kink" energy of field configurations. There- 
fore the kinetic energy in Eq.(69) suppresses the 
spatial dependence of the field configuration B^J^^ 
for small i?, i.e., for high density state, giving flat 
configuration in the high density phase. 

2) The mass term m^tr{i3}j"''i3}j"''} suppresses the ab- 

solute value of field more severely in accor- 
dance with its larger mass m^. 

3) The couplings between pions and heavier (axial) 
vector mesons i?^"' with larger index n are found to 
become smaller, which is suggested by the "bottom- 
up" dimensional deconstruction model [31] and also 
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FIG. 2: Energy density of single baryon for the Skyrme model 
with radius i? of 5*^. Below the critical radius, R < R^^^^"^" = 
V^, energy density becomes a uniform distribution as the 
identity map, shown by the dashed lines. 



FIG. 3: Energy density of single baryon for the BIS model 
with radius ii of S^. Below the critical radius, R < Rf^ = 
1.19, energy density becomes a uniform distribution as the 
identity map, shown by the dashed lines. 



the "top-down" holographic approach [9]. In fact, 
within the holographic models, the origin of the 
meson mass in four dimensional space-time can be 
regarded as the oscillation of meson wave function 
in the extra fifth dimension z, denoted by the mass 
relation nvf^ = A„ in Sec. II. Such larger oscil- 
lations of heavier (axial) vector mesons wave func- 
tions in the fifth dimension have the smaller overlap 
with that of pions, giving the smaller coupling con- 
stants with pions in four dimensional space-time [9] . 
In fact, some recent experiments with the hadron 
reactions provide some interesting data, showing 
that the heavier (axial) vector mesons tend to have 
smaller width for the decay into pions despite of 
larger phase space [52], which may be consistent 
with the prediction of holographic QCD as their 
smaller coupling constants with pions as mentioned 
above. Hence the effects of heavier (axial) vector 
mesons should be smaller for a baryon large- 
amplitude pion field, i.e., the chiral soliton. 

These considerations 1), 2), 3) about the meson effective 
action in holographic QCD should support our conjecture 
about the general disappearance of (axial) vector meson 
fields in the high density phase. In other words, only the 
pion fields survive near the critical density in the large- 
Nc baryonic matter. We call this phenomena as "pion 
dominance" near the critical density. 

In Sec. VI, we will show that, even if the p meson field 
disappears near the critical point, it affects the critical 
density of the phase transition through its mass and also 
the interactions with the pion field in the action S'ofr in 
Eq.(19). In this sense, even with the pion dominance 
near the critical point proposed above, the (axial) vector 
mesons may still affect the critical phenomena like the 
critical density through their contributions in the effec- 
tive action if they arc included. 



As for the survival of the pion fields near the criti- 
cal density as the pion dominance, we give the follow- 
ing reason: the unit baryon-number on the manifold 
comes from the boundary conditions of the pion hedge- 
hog configuration F{r) at the north (r = 0) and the 
south (r = ttR) poles on as F{0) = tt and F(7ri?) = 
in (62). In this sense, pions play the essential roles for 
the baryon-number current with their field boundaries. 
Therefore the pion field cannot disappear because of the 
baryon-number constraint (62) for each unit cell of the 
manifold S^. On the other hand, there is no constraints 
for the other (axial) vector meson fields, and they can dis- 
appear in the high density phase as some representation 
of "deconfinement" . 

Such dominance of the pion field near the transition 
point might be somehow related with the appearance of 
chiral plasma modes tt and a above Tc of the decon- 
finement phase transition observed in the lattice QCD 
stuty [53]. In fact, the screening masses were measured 
for a variety of color-singlet channels with the quantum 
number of tt, cr, p, 6i, and ai mesons. The chiral multi- 
plet of TT and a was found to appear as the bound state 
even above whereas the other (axial) vector mesons 
would be in the continuum with the threshold as the twice 
of the lowest Matsubara frequency, i.e., 27rr. Such sur- 
vival of the chiral plasma modes above Tc would inspire 
the concept of the strong-coupling quark-gluon plasma 
(sQGP) due to the long-range, non-perturbative effects 
of the strong interaction even at the highest temperature, 
instead of the naive free gas picture of quarks and gluons 
only due to the asymptotic freedom. Such strong corre- 
lation would also develop even at dense regime of QCD, 
which may give some linking of our analysis with "pion 
dominance" in dense QCD. 

By seeing some QCD phenomenologies, one can also 
find that quark degrees of freedom arc often represented 
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FIG. 4: Energy density of single baryon and p meson contri- 
butions in the terms of the BIS model (64) for R = 4.0, \/2, 
and 1.19. Labels for lines in (a) denote each term in (64), e. 
g. "4p" corresponds to (— l/2)(;4p / d*xtr[p^, p^]^ (See Ref [9] 
for labels). The vertical dashed lines in (b) and (c) show the 
boundaries of at the south pole r = nR, respectively. At 
the critical radius, 7? = -R^'f = 1.19 in (c), energy density be- 
comes a uniform distribution as the identity map, where the 
p meson contributions disappear and only pion contribution 
survives as the "pion dominance". 



by "pions" . For example, in the case of the chiral quark 
model [54, 55], quark-antiquark correlations are repre- 
sented by pionic collective modes through the bosoniza- 
tion scheme for QCD, giving the "Cheshire cat picture" 
where the quark dynamics are represented by pions. In 
fact, it have provided a foundation of Skyrmc soliton pic- 
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FIG. 5: Pion field F{r) of the BIS model for R = 4.0, ^/2, and 
1.19. At the critical radius, R = Rf^^ = 1.19, F{r) coincides 
with a linear configuration as the identity map in Eq. (67). 



(Rescaled p-meson field G(r) for 7?=4.a2'^, 1. 19) 
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FIG. 6: Rescaled p meson field G{r) = ^ of the BIS model 

for R - 4.0, V2, and 1.19. At the critical radius, R = Rf,^^ = 
1.19, G(r) becomes zero configuration as the identity map in 
Eq. (67). 



ture for the baryon as a large-amplitude pion field. These 
traditional phenomenologies may also support the domi- 
nance of the pions near the critical density in the baryonic 
matter, with the appearance of quark-gluon dynamics. 



C. Mass and root-mean-square mass radius 

In this section, wc discuss the baryon-numbcr density 
dependence of the mass and the root-mcan-square (RMS) 
mass radius of single baryon for the Skyrme model and 
the BIS model, by changing the size of the manifold S^. 

In Fig. 7, we show the mass of the baryon for the 
Skyrme model and the BIS model with radius R oi S^. 
In the case of the Skyrme model, the mass decreases from 
the value on the flat space : 1.231 x 12tt^[^] [45] as 
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(Mass of a baryon with radius R) (RMS mass radius with radius R) 




FIG. 7: Mass of single baryon for the identity map, the 
Skyrme model, and the BIS model with radius R o{ S^. The 
variables "1.231" and "1.115" are the masses of the baryon on 
the flat space R'' for the Skyrme model and the BIS model, 
respectively. 



FIG. 8: RMS mass radius of single baryon for the identity 
map, the Skyrme model, and the BIS model with radius R 
of S^. The variables "1.422" and "1.268" are the RMS mass 
radius of the baryon on the flat space for the Skyrme 
model and the BIS model, respectively. 



radius R as 



R decreases, and it coincides with the mass of identity 
map below the critical radius, R < R^^{["^'^ = \/2 [27]. 
This coincidence of numerical solution with identity map 
indicates the transition from the localized phase to the 
uniform phase as the delocalization phase transition. In 
the case of the BIS model, the mass decreases from the 
value on R'^ : 1.115 x 127r^[||] [9] as R decreases, and it 
coincides with the mass of identity map below the critical 
radius, R < R^^^ = 1.19, which is smaller than 7?^^!^^'"^'' 
due to the shrinkage of baryon by the p meson effects. 

From a physical point of view, the decrease of the 
baryon mass as R decreases in Fig. 7 might come from 
the partial restoration of chiral symmetry in the dense 
baryonic matter, which will be discussed with the order 
parameter in Sec. V. This result is related with the well- 
known picture of the baryon mass generation from the 
quark condensate in the vacuum, suggested in some QCD 
phenomenologies, e.g., the QCD sum rules with loffe for- 
mula [56] , and also other chiral effective models [57, 58] . 
In fact, the decrease of the baryon mass with the chiral 
symmetry restoration is observed in the finite tempera- 
ture lattice QCD study [42]. Note also that the decrease 
of the baryon mass with the chiral symmetry restoration 
in dense QCD is proposed in the framework of AdS/QCD 
model as the phcnomcnological bottom-up constructions, 
called "hard- wall" model [37]. 

Next we analyze of the RMS mass radius of the baryon 
for the Skyrme model and the BIS model. By using the 
normalized energy density e(r) = e{r)/E {eijr) is total 
energy density and E is the mass of single baryon), the 
RMS mass radius can be naturally introduced on with 



-,1/2 



r 

I iirdrR^ sin^ — • e(r)r^ 
Jo R 



1/2 



, (70) 



where R^ sit? denotes the measure of the curved man- 
ifold . In the case of the identity map, the RMS mass 
radius can be analytically calculated by using the nor- 
malized uniform energy density for the identity map : 
ffid = l/27r2i?3 as 



AndrR^ sir? — ■ Sidr'^ 

R 



1/2 



(71) 



In Fig. 8, we show the RMS mass radius for the Skyrme 
model and the BIS model with radius i? of S''^. One can 
find that the RMS mass radius in each model increases 
non-linearly around its critical radius as R decreases, in- 
dicating the "swelling" phenomenon of the baryon in high 
density phase. Such swelling of the baryon can also be 
seen in the analysis of the Skyrme crystal on the three- 
dimensional cubic lattice [26] and also the finite density 
bag model [59]. After the delocalization phase transition, 
the energy density of single baryon is uniformly saturated 
over the surface of S^. Then the RMS mass radius de- 
creases linearly with the decrease of R as seen in Eq.(71). 
(In the uniform phase, baryons are no longer localized, 
and the RMS mass radius may have less physical mean- 
ing.) 
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D. Swelling mechanism and its phenomenological 
implications 

In the previous section, we numerically find the 
swelling of the baryon by taking the Skyrmion picture for 
the baryon analysis. Actually, as discussed in Sec. Ill A, 
the Skyrmion picture is naturally derived from QCD 
through the new concept of holography. Therefore one 
can expect the swelling phenomenon in QCD itself. In 
this section, we try to give a qualitative explanation 
about the swelling mechanism of a baryon in terms of 
QCD, by especially referring to contexts of the lattice 
QCD [42] and also the finite density bag picture [59]. 

First, in lattice QCD study, the non-perturbative QCD 
vacuum is analyzed within a finite box of periodic bound- 
aries (and also Dirichlct boundaries) in the Euclidean 
space. The non-perturbative effects with long wave 
length tends to be "blocked out" from such finite box 
of definite boundaries, so that the QCD vacuum in the 
finite box is found to approach to the perturbative one. 
Also, the finite-temperature phase transition in lattice 
QCD occurs because of the definite boundaries along 
the imaginary time axis in Euclidean space; as the tem- 
perature increases, the QCD vacuum in the Euclidean 
space is more closely packed along the imaginary time 
axis to block out the non-perturbative effects with long 
wave length, which eventually gives the transition into 
the perturbative QCD vacuum, like the chiral symmetry 
restoration and also the deconfinement. 

Next we combine the above picture into the finite den- 
sity bag model, to explain the mechanism of the swelling 
in QCD. In the bag model, the baryon is represented 
by a "bag" with the perturbative QCD vacuum sur- 
rounded by the Dirichlet boundaries, which is supported 
by the non-perturbative QCD vacuum with the bag pres- 
sure as shown in Fig. 9(a). As a number of bags in- 
creases, representing the high density baryonic matter, 
the non-perturbative QCD vacuum around the bags is 
more closely packed within a definite region, imitated 
by a cube with dashed lines in Fig. 9(b). Then, in 
analogy with the case of lattice QCD study discussed 
above, the QCD vacuum around the bags should ap- 
proach to the perturbative one by blocking out the non- 
perturbative effects with long wave length as shown by 
the several waves in Fig. 9. Therefore the bag pressure 
as a non-perturbative effect decreases as the number of 
bags increases, so that it eventually gives the swelling 
of bags. Such swelling in turn decreases the region of 
non-perturbative QCD vacuum around the bags to give 
the decrease of the bag pressure again. In this sense, 
such swelling could occur non-linearly as shown in Fig. 8, 
giving some drastic change of features like "deconfine- 
ment" in the baryonic matter. One can expect some 
resemblance between the finite density bag model and 
the Skyrme matter in our study with respect to bary- 
onic matter, so that the swelling phenomena as shown in 
Fig. 8 could be understood on the same footing provided 
above. 



Now we also suggest the physical effects of the swelling 
in the high density baryonic matter. Here we take into 
account possible baryon excitation with the chiral soliton 
picture. Since a moment of inertia of the baryon is related 
with the size of the baryon, some non-linear increase in 
the moment of inertia around critical density is expected 
along with the swelling of a baryon. In the semiclassical 
quantization procedure of a Skyrme soliton, the baryon 
mass spectra is given [43] as 
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(72) 



where Mhh is static hedgehog mass and T is the mo- 
ment of inertia of a baryon. The quantum number J 
in Eq. (72) denotes the spin S and the isospin / of the 
baryon as J{= S = I) = ^ for N and J{= S = I) = ^ 
for A. The baryon mass splitting in the second term in 
Eq.(72) is proportional to the inverse of the moment of 
inertia, indicating that a larger object cannot be easily 
rotated. Therefore we can now propose that the baryon 
mass splitting decreases as the density increases in the 
medium, because the moment of inertia X increases due 
to the swelling. In other words, the individuality of each 
baryon like TV, A, etc., would be lost at least with respect 
to their mass spectra as the density increases. Such lost 
of individuality may be regarded as some precursory phe- 
nomena of deconfinement, where the dominant degrees of 
freedom in the hadronic matter shifts from the confined 
composites as baryons (and mesons) into quarks and glu- 
ons. 

Such swelling of the baryon leads to some interest- 
ing phenomenological realizations. As one possibility, 
here we propose the stable N-A mixed matter in dense 
QCD. For simplicity, we consider the "symmetric nu- 
clear matter" , where only the strong interaction is taken 
into account without the electro-magnetic interaction. In 
the low density case, there exists large A^-A splitting as 
{Ma -Mn) ^ 290MeV, and the Fermi surface of the nu- 
cleons fiN may locate below the threshold of A isobars, 
i.e., its mass Ma. This situation gives only the nucleon 
degrees of freedom as nuclear matter. Now, if the swelling 
occurs as the density increases, the N-A splitting should 
decrease as discussed in Eq. (72). Then, the Fermi sur- 
face of the nucleons /xjv would excess Ma to give the 
A isobar degrees of freedom in the medium through the 
equilibrium process A^ -I- tt ^ A. In this sense, A^-A 
mixed matter would be realized in dense QCD because 
of the swelling of a baryon. Such A^-A mixed matter 
may appear in the deep interior of neutron stars between 
the nuclear crust and the core of quark matter as the 
precursor of deconfinement, giving some softening of the 
EOS of neutron stars relative to the analysis without the 
mixed matter. While these are qualitative discussions, 
we feel that these explanations catch some essential as- 
pects of baryonic matter with baryon excitation in view 
of its "swelling" nature. 
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FIG. 9: Schematic figure of the finite density bag model for low density case (a) and high density case (b): the baryon is 
represented by a "bag" with perturbative QCD vacuum (P-QCD) surrounded by Dirichlet boundaries, supported by the non- 
perturbative QCD vacuum (NP-QCD) with the bag pressure B. Several waves in (a) and (b) represent the non-perturbative 
effects of QCD with certain wave length. A cube with dashed lines in (b) imitates a definite region of non-perturbative QCD 
vacuum surrounded by the surfaces of the bags with Dirichlet boundaries, to discuss the mechanism of "swelling" in the high 
density baryonic matter. 



V. PHASE TRANSITIONS WITH ORDER 
PARAMETERS 

In this section we discuss the delocalization phase tran- 
sition in view of the deconfinement, and also the chiral 
symmetry restoration by introducing the proper order 
parameters in the Skyrme model and the BIS model. The 
relations between these phase transitions are also consid- 
ered by referring to other QCD phenomenologies in the 
end of Sec. VB. 



A. Delocalization phase transition 

First we consider the delocalization phase transition 
with the order parameter. By using the normalized en- 
ergy density e(x) = e{K)/E (e(x) is total energy density 
and E is the mass of single Skyrmion), one can introduce 
the spatial fluctuation of the energy density around 
the uniform energy density distribution [60] as 



--id\ 



(73) 



where Sid 



is the normalized uniform energy den- 



sity for the identity map. If the Skyrmion is well localized 
like the delta function, the energy density e(x) is almost 
decoupled with uniform energy density distribution £id 
through the spatial integral over 5''^. Therefore its spa- 
tial fluctuation $(i?) becomes. 



<f>(i?) ^i/53rf3x{|e-(x)| + |e-id|} = l. 



(74) 



On the other hand, if the Skyrmion gives uniform en- 
ergy density distribution as the identity map, the energy 
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FIG. 10: Spatial fluctuation of the energy density of 

single baryon for the Skyrme model and the BIS model with 
radius _R of S^. 



density £(x) coincides with £id, giving 
$(i?) = 0. 



(75) 



In this sense, $(i?) in Eq.(73) measures the amount 
of localization in the energy density distribution of the 
Skyrmion, regarded as the order parameter of delocaliza- 
tion phase transition in the baryonic matter. 

In Fig. 10, we show the value of ^{R) for the Skyrme 
model and the BIS model with radius R of S*^. For suf- 
ficiently large R as the low density state, the Skyrme 
soliton is well localized for both the Skyrme model and 
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the BIS model; giving $(i?) ~ 1. With the decrease 
of R, the baryon tends to delocahze with the decrease 
of ^{R), regarded as some medium effect in the bary- 
onic matter. The delocahzation phase transition into 
the uniform phase with <i>(i?) = occurs at the criti- 
cal radius R = R^^^'^™'^ = \f2 for the Skyrmc model and 
at i? = R^f^ =1.19 for the BIS model. Actually, the 
delocahzation phase transition in the BIS model is de- 
layed along with the decrease of R because of the heavy 
p mesons in the core region of the baryon, discussed in 
Sec. IV A. 

Here one may have some temptation to relate the de- 
localization phase transition with the "deconfinement" 
of the baryon in QCD, which should be carefully dis- 
cussed. Actually, there have been several conflicts about 
the meaning of the delocahzation phase transition in the 
standard Skyrme model with finite density [26, 27, 60- 
62]. Within the traditional works in 1970', the direct 
linking of the Skyrme model with QCD was still un- 
certain. Furthermore, the standard Skyrme model does 
not manifestly include quarks and gluons, so that large 
care should be taken to relate the delocahzation of the 
Skyrmion with the deconfinement of the baryon as the 
appearance of quark-gluon dynamics. In our case, how- 
ever, the new concept of holography allows us to derive 
the Skyrmion picture from QCD itself, which is discussed 
in Sec. Ill A. Moreover, one should also note that quark 
degrees of freedom can be represented by "pions" , sug- 
gested in some QCD phenomenologies with "Cheshire cat 
picture" as in the chiral quark model [54, 55] derived 
from the bosonization scheme for QCD. With these back- 
grounds, here we propose that the delocahzation phase 
transition in the BIS model can be more admissibly re- 
lated with the deconfinement of the baryons in QCD, 
relative to the standard Skyrme models in 1970'. 



B. Chiral symmetry restoration in non-linear 
realization 



Next we consider the chiral symmetry restoration with 
the order parameter. Normally, in case of the linear 
sigma model, the chiral symmetry restoration can be sig- 
naled by the vanishing of the expectation value of the 
meson fields at the matter ground state as 



are introduced from the chiral field U{x^) G SU(A^/)a as 



(76) 



where the sigma meson field cr(a;^) and three pseu- 
doscalar pion fields n(a;^) are called as the "linear" real- 
ization of chiral symmetry. We relate this terminology to 
the non-linear sigma model, since the meson effective ac- 
tion and also the Skyrmion picture are based on the "non- 
linear" realization of the chiral symmetry. In the non- 
linear sigma model, the meson fields cr(x^) and n(x^) 



= cos{|7r|//^} + iTaTTa sin{|7r|//^} 

= {<T{Xf,) +iT -UiXf,)}//^. 

{\n\= ^^Tra7ra,Tra = TTa/\Tr\) (77) 

Since U{Xf^) is unitary matrix with the condition 
U'^{x^)U{x^) = 1, the squared sum of the meson fields 
<j{x^) and n(a;^) in Eq. (77) should be constant every- 
where as 



(78) 



Therefore, as far as the action is written in the chiral field 
U{x^) with fixed /,r, the meson fields cr(a;^) and Tl{x^) 
arc forced on a surface of a three-dimensional closed man- 
ifold S'fjjj with finite radius in the internal space. In 
fact, the existence of such closed manifold S^^^ is essen- 
tial for the concept of a Skyrmion, which is a non-trivial 
winding of the compactified physical manifold around 
the other closed manifold 5"?^^ with conserved topolog- 
ical charge, belonging to the homotopical classification 
"■3(5''^) ~ Z [63]. However, by comparing Eq. (78) with 
Eq. (76), one would always encounter a problem of how 
to describe the chiral symmetry restoration with its non- 
linear realization. 

Now in this paper, we take a spatially-averaged con- 
densate of the meson fields cr(x^) and n(a:^) as the order 
parameter of the chiral symmetry restoration with the 
non-linear realization [62] as 



{{a{x,)Y + {IV{:x,,)Y}/fl. 



(79) 



Here, the bracket in Eq.(79) denotes the three- 
dimensional spatial average in the medium with volume 
V as 



(©(%)> ^ 
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(80) 



Such spatially-averaged condensate of the meson fields 
is somehow similar to the spatially-averaged magnetiza- 
tion (M(x)) as the global order parameter of the ferro- 
magnetic material like bulk iron. Within the bulk iron 
at sufficiently high temperature, the spins of the ions 
orient randomly to get the entropy gain in the free en- 
ergy. As the temperature decreases below the critical 
temperature without the external magnetic field, there 
appear magnetic domains with non-zero local magneti- 
zation M(x) 7^ 0. However, because of the total angular 
momentum conservation, the spatially-averaged magne- 
tization (M(x)) should vanish with the appearance of a 
complex structure of magnetic domain walls within the 
bulk iron, which gives no breaking of the global spa- 
tial symmetry macroscopically in the bulk material as 
(M(x)) = 0. In this sense, the spatially-averaged con- 
densate of the meson fields in Eq.(79) can be regarded as 
the "global" order parameter of the chiral symmetry in 
the bulk hadronic matter. 
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By taking the hedgehog configuration Ansatz 
T^aiXfj^)/ /tt = XaFir) as in Eq.(35), the meson fields 
(y{x^) and n(x^) can be written as 

<yix^)/U = cos{7r(a;^)//^} = cos i^(r), (81) S 

n(x,J//^ = 7rsin{7r(x^)//^} = xsinF(r). (82) o 

a 

Note here that the pion field with the hedgehog ansatz in ^ 
Eq.(82) is proportional to the unit directional vector x, ^ 
so that its spatial average becomes trivially zero. There- ^ 
fore only the sigma meson field should be considered as 
the global order parameter of the chiral symmetry. If a 
Skyrmion is well localized around the north pole at r = 
on the manifold S^, the classical meson filed configura- 
tion F(r) becomes zero almost everywhere except for the 
localized point r ^ 0. Therefore the sigma meson config- 
uration in Eq.(81) becomes aix^^j ^-^ ~ 1 almost every- 
where. In other words, the vacuum is so much oriented 
even globally in the medium, thus its spatial average be- 
comes non-zero as 
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FIG. 11: Spatially-averaged condensate of the sigma meson 
field {(^{x^)) / fjv of single baryon for the Skyrme model and 
the BIS model with radius _R of S^. 



(a(.T^))//. ~ 1, 



(83) 



indicating the appearance of the chiral symmetry bro- 
ken phase. On the other hand, if the Skyrmion gives 
uniform energy density distribution as the identity map 
F{r) = TT — r/R in Eq.(67), the sigma meson configura- 
tion becomes a{Xf^)/ f^r = cos changing monotonously 
from —1 to 1 over the coordinate space through the 
arc distance r G [0,7ri?]. Therefore its spatial average 
vanishes as 



(fT(x^))//. = 0, 



(84) 



indicating the appearance of the chiral symmetry re- 
stored phase. These considerations about Eqs. (83) and 
(84) also indicate that the chiral symmetry restoration 
is indirectly related with the energy density distribution 
in the hadronic matter through the classical meson field 
configurations. 

In Fig. 11, we show the spatially-averaged condensate 
{(j{x^)) If-,,: for the Skyrme model and the BIS model with 
radius i? of 5''^. The spatial average with the hedgehog 
configuration (81) is explicitly taken for each closed man- 
ifold as 



= {cos Fir)) 

/g3 a; COS i^(r) 

J^" AirdrR^ siv? • cosF(r) 
27r2i?3 ■ 



(85) 



For sufficiently large R as the low density state, the 
Skyrme soliton is well localized for both the Skyrme 
model and the BIS model, as previously shown in Fig 10. 
Therefore, as discussed in Eq. (83), the spatially-averaged 
condensate of the meson field becomes (cr(x^t)) //^r ~ 1 as 
the chiral symmetry breaking in both models as shown 



in Fig. 11. In this sense, p meson field (and also the 
other vector and axial vector meson fields, if they are 
included) would have less importance for the chiral sym- 
metry breaking in the bulk hadronic matter as far as suf- 
ficiently localized baryon appears in the medium. With 
the decrease of R, (cr(x^))//7r monotonously decreases, 
which corresponds to the partial chiral symmetry restora- 
tion as the medium effect in the baryonic matter. The 
chiral symmetry is fully restored with (17(0;^))//^ = at 
the critical radius R ~ — for the Skyrme 

model and at i? = R^^^^ =1.19 for the BIS model. Ac- 
tually the small radius R, i.e., the larger baryon- number 
density is needed to give the chiral symmetry restoration 
due to the p meson fields, in similar to the case of the 
delocalization phase transition as previously discussed. 

Finally, we comment about the relation between the 
deconfinement and the chiral symmetry restoration in 
QCD, which has been discussed so long in these decades 
of the hadron physics. The relation between two phase 
transitions is nontrivial since they are characterized by 
the different symmetries: the global chiral symmetry 
SU(A^/)l X SU(A^/)r, and Z(3) symmetry as the cen- 
ter of SU(3)c gauge group. Note that the former can be 
defined in the massless quark limit, while the latter in 
the heavy quark limit, leading to the difficulty to discuss 
their relation. 

Despite of the difficulty, several implications have been 
obtained: In finite temperature, the lattice QCD studies 
suggest that the two phase transitions occur at the same 
critical temperature Tc ~ 170MeV [42]. Such simul- 
taneous occurrence is supported by the analyses of the 
Nambu- Jona-Lasinio model with a Polyakov loop (PNJL 
model) as a low-energy effective theory of QCD [64]. Ac- 
tually, these two phase transitions are separated with re- 
spect to their dominant symmetries as mentioned above, 
so that such mysterious coherence between two phase 
transitions is often compared to some kinds of "entangle- 
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mcnt" [64]. Furthermore, the simultaneous occurrence of 
the two phase transitions is also suggested by the recent 
analysis of the holographic QCD with D4/D8/D8 multi- 
D brane configurations: two independent chirality spaces 
on the D8 and D8 branes respectively are connected with 
each other (chiral symmetry breaking) by the "worm 
hole" of the D4 supergravity background into which the 
colored objects like quarks and gluons arc absorbed (color 
confinement) [9]. In this sense, these two phase transi- 
tions are more directly related with each other as just sin- 
gle event on the "worm hole" in the extra-dimensions, the 
effects of which might appear as some kinds of mysterious 
entanglement in view of the four-dimensional space-time. 

Now, the situation has not been clear in finite density. 
In this work, we indirectly relate the two phase transi- 
tions through the meson field configurations. Admittedly 
our approach includes only meson fields and baryons ap- 
pear as mesonic solitons. However, if the Cheshire cat 
picture holds and quark-gluon dynamics can be indi- 
rectly expressed in the meson dynamics as an effective 
model, our results could be interpreted as the simultane- 
ous occurrence of two phase transition in the finite den- 
sity QCD. Investigations in this direction arc interesting 
and should be developed further. 
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FIG. 12: Total baryon-number density pB in Eq. (86) with 
radius i? of S"^, where the experimental inputs are taken as 
= 92.4MeV and e = 7.315 in Eqs. (87) and (88). The ver- 
tical axis is given as pB divided by the normal nuclear density 
po — 0.17fm~''. Critical densities of the phase transitions for 
the Skyrme model and the BIS model are also shown for each 
critical radius as -R^^f = \/2 and = 1.19. 



VI. CRITICAL DENSITY WITH PHYSICAL 
UNITS 

In this section, we show the critical densities of the 
phase transitions for the Skyrmc model and the BIS 
model with recovering the physical units. 

In this study, single baryon is placed on a closed mani- 
fold with the surface volume 27t^R^, so that the total 
baryon-number density ps can be given as 

ps = (27r2i?3)-i(e/^)3[fi^-3]. (86) 

Now, in the holographic model, by fixing two parame- 
ters, e.g., experimental inputs for and irip, all the phys- 
ical quantities like masses and the coupling constants are 
uniquely determined, which is a remarkable consequence 
of the holographic approach discussed in Sec. II. First, 
we take and mp as experimental values, 

U = 92.4MeV, rup = 776.0MeV. (87) 

Then the Skyrme parameter e in Eq.(42) can be uniquely 
determined as 

e~ 7.315. (88) 

With these experimental inputs, we show in Fig. 12 the 
total baryon-number density in Eq- (86) as a function 
of radius i? of S*^, divided by the normal nuclear density 
Po — 0.17fm~^. From Fig. 12 and also Eq.(86), the critical 
densities for the Skyrme model (i?^jl^j^""° = v^) and the 



BIS model (i?^rU = 1-19) ^''^ found as follows: 

Pb''"" ^ {stt^ (^-r")'} \e/.)^ 4.26po(89) 

pT ^ {27r' {R^Mf}'' ^'M' =^ ^-l^Po. (90) 

The heavy p meson fields in the core region of the baryon 
tend to decrease the total size of the baryon [9], so that 
the p meson field has a significant role to increase the 
critical density into the uniform phase as in Eq. (90) rel- 
ative to Eq. (89). 

In Sec. IV A, we discuss the disappearance of the p me- 
son field in high density phase of baryonic matter, while 
we now find the significant roles of p meson field as for 
the critical density in Eq. (90) relative to Eq. (89). This 
situation is somehow resemble to the "two-Higgs model" 
with scalar fields (j) and x in the finite temperature. On 
the critical temperature of a phase transition with the 
dynamics of a scalar filed X: the condensate of the x held 
might become trivial as (x) = 0. However, in general, the 
mass of X field and also the interactions between x and 4> 
fields affect the critical phenomena like critical tempera- 
ture. In this sense, even if p meson field disappear near 
critical point, they could affect the critical phenomena 
as in Eq. (90) through its mass and also the interactions 
with pions in the action Scs in Eq. (19). 

Here we briefiy comment about the value of the critical 
density by taking other parameter set used by Adkins 
et al. [25]. In 1983, Adkins et al. analyzed the baryon 
mass spectra with the semiclassical quantization of the 
Skyrmc model, whereas the A^-A splitting corresponds 
to the higher order contribution 0{N~^) relative to the 
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static Skyrme action O(iV^) in the large- A^c expansion. 
To reproduce the proper amount of the N-A sphtting, 
they take the smaller pion decay constant and Skyrme 
parameter [25] as 

U ^ 64.5MeV, e = 5.44. (91) 

Now if the parameter set in Eq.(91) is used for the 
baryon- number density in Eq. (86), we find too small crit- 
ical density (^ 0.6po) for the Skyrme model. If one can 
relate the phase transition in the Skyrme model with the 
deconfinement of the baryon in QCD, it might propose, 
e.g., that quark degrees of freedom are already mani- 
fest in the normal nuclei. This unreasonable results may 
come from the fact that scmiclassical quantization proce- 
dure corresponds to extracting higher order contributions 
0{N~^) from the leading order Skyrme action 0{N^) as 
for the large- TVc expansion. 

Same comments can also be applied to the holographic 
model if one try to perform the scmiclassical quantiza- 
tion as the baryon analysis. In fact, holographic QCD is 
derived as the large- iVc effective theory even by starting 
from superstring framework. Higher order effects of the 
large-TVc expansion like the baryon mass splitting cor- 
respond to the string loop effects beyond the classical 
supergravity, which should be fairly untractable. With 
these considerations above, we do not intentionally pro- 
ceed to the scmiclassical quantization of Skyrme soliton 
in the present work, and we employ the experimental in- 
puts (87) and Skyrme parameter (88) in the analysis of 
the critical densities with recovering the physical units. 
Actually, due to the scaling property of the BIS model 
discussed in Sec. II, all the results in ANW unit in the 
previous sections should not be altered, being indepen- 
dent of the definite values for and nip. 

VII. SUMMARY AND OUTLOOK 

We have studied baryonic matter in holographic QCD 
with D4/D8/D8 multi-D brane configurations, by analyz- 
ing the system of single branc-induced Skyrmion on the 
three-dimensional closed manifold S^. By changing the 
size of S^, the density dependence of the baryon proper- 
ties are examined from a holographic point of view. 

First we begin with the Dirac-Born-Infcld (DBI) ac- 
tion of the probe D8 brane with D4 supergravity back- 
ground. With the dimensional reductions, we get the 
five-dimensional Yang-Mills action with a curved fifth di- 
mension, as the leading order of the large-A^c and large 
'tHooft coupling expansions in dual of non-pcrturbative 
(strong-coupling) QCD. Through the mode expansions 
of the five-dimensional gauge fields, we get the four- 
dimensional meson effective action from holographic 
QCD. Especially, we emphasize the appearance of the 
ultraviolet cutoff scale A/kk ^ IGeV in the holographic 
approach to be dual of QCD. Therefore we construct the 
four-dimensional meson effective action with pion and p 
meson fields below the cutoff scale Mkk- 



Then we discuss the baryon in holographic QCD as 
the "brane-induced Skyrmion" in the four-dimensional 
meson effective action. The analyses of the baryon in 
the meson effective action with restricted degrees of free- 
dom below the cutoff scale Mkk is called the "truncated 
resonance model" for the baryon. Taking the hedgehog 
configuration Ansatz for pion and p meson fields, we can 
investigate many properties of single baryon as the brane- 
induced Skyrmion, which arc inclusively summarized in 
Ref. [9]. 

Then we consider the baryonic matter in holographic 
QCD, as the extension of holographic approach to dense 
QCD. Especially we treat the baryonic matter with large- 
Nc, because holographic QCD is derived as the large- A^c 
effective theory, in dual of the classical supergravity. For 
sufficiently large Nc, the kinetic energy and the quantum 
effects can be suppressed relative to the static mass from 
simple large- A'^c countings [8], where the baryonic mat- 
ter comes into the static Skyrme matter. In this sense, 
we analyze the static Skyrme matter to see the typical 
features of the baryonic matter with large- TVc conditions. 

In order to analyze the static Skyrme matter on the 
flat coordinate space R^, we alternately treat the sys- 
tem of single brane-induced Skyrmion on the three- 
dimensional closed manifold S^. The interactions be- 
tween the baryons arc simulated by the curvature of the 
closed manifold S^, and, as the size of decreases, the 
baryon-number density increases in this modeling. Ac- 
tually, through the projection procedure from the flat 
space R"^ onto the curved space S^, we get the hedgehog 
mass and also the Euler-Lagrange equations for pion and 
p meson fields of the brane-induced Skyrmion on S^. 

By numerically solving the Euler-Lagrangc equations 
for pion and p meson fields on S^, we find a stable soliton 
solution as the brane-induced Skyrmion on . By using 
this solution, we analyze many properties of the baryon 
within the baryonic matter. Especially by comparing the 
standard Skyrme model without p mesons and the brane- 
induced Skyrme (BIS) model, the roles of (axial) vector 
mesons in dense QCD are discussed from a holographic 
point of view. 

First we show the baryon-number density dependence 
of the energy density and the field configuration profiles 
of single baryon, by changing the size of S^. As the size 
of decreases, the localized energy density distribu- 
tion of single baryon becomes uniform one as the iden- 
tity map, which is called the "delocalization phase tran- 
sition" . The critical radii of the phase transitions are 
given as i? = R^m'""" = for the Skyrme model 

and R = Rf^ = 1.19[^] for the BIS model. Because of 
the shrinkage of the baryon size due to the p meson ef- 
fects [9], the smaller critical radius of S^, i.e., the larger 
baryon-number density is needed for the BIS model to 
give the delocalization phase transition. We also find that 
p meson field absolutely disappears and only pion field 
survives near the critical density. Then, with the math- 
ematical arguments, we propose a remarkable conjecture 
that all the (axial) vector meson fields would disappear 
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and only pion field survives near the critical density, re- 
ferred as the "pion dominance" in dense baryonic matter. 

We also investigated the baryon-number density de- 
pendence of the mass and the root-mean-square (RMS) 
mass radius of single baryon. We find the the "swelling" 
phenomenon of the baryon, as the non-linear increase of 
the RMS mass radius near the critical density. Actually, 
non-perturbative QCD vacuum around the baryons are 
closely packed as the baryon-number density increases. 
Therefore, the non-perturbative effect with long wave 
length like the bag pressure would be blocked out from 
such definite region to give the swelling of the baryon. 
Such swelling provides the the decrease of the baryon 
mass splitting, indicating the lost of the individualities 
for the baryons as the precursor of the deconfinement. 
We also propose the stable N-A mixed matter in dense 
QCD because of the swelling, which could, e.g., soften 
the EOS of the neutron stars. 

The features of the delocalization phase transition and 
the chiral symmetry restoration are also analyzed with 
the order parameters in the holographic QCD. We conjec- 
ture with careful arguments that the delocalization phase 
transition can be related with the deconfinement of the 
baryons in QCD. Then we find the coherence of the de- 
confinement and the chiral symmetry restoration through 
the meson field configurations. Such coherence of the two 
phase transitions are also suggested in the lattice QCD 
study in the finite temperature [42], the PNJL model as 
a low-energy effective theory of QCD [64], and also the 
Sakai-Sugimoto model as one of the reliable holographic 
approaches [9]. 

We also calculate the critical densities of the phase 
transitions with the experimental inputs for the pion de- 
cay constant and p meson mass rup as = 92.4MeV 
and nip = 776.0MeV. We find the critical densities as 
^^kyrmc ^ 4.26/90 for the Skyrme model and p^^^ ~ 
7.12/90 for the BIS model. We can see that the larger 
baryon-number density is needed for the BIS model to 
give the phase transitions, because of the shrinkage of 
total size of the baryon with heavy p mesons in its core 
region. 

Finally we compare our truncated resonance approach 
for the baryons with the other works of baryons as the 
instantons in holographic QCD. In fact, there seem to 
exist some conflicts for the baryon analysis in holo- 
graphic QCD, especially between the truncated reso- 
nance model and the instanton models. In Refs. [11-14], 
the baryons are studied as the instantons on the five- 
dimensional gauge theory of the probe D8 brane with 
D4 supergravity background as the holographic dual of 
QCD. The properties of the baryonic matter are also 
analyzed by the system of the single instanton on the 
three-dimensional closed manifold as the Wigner-Seiz 
approximation [15]. In these analyses, the instanton is 
found to shrink into zero size only for the DBI sector of 
the effective action of the D8 brane, which corresponds to 
the leading order of large- A'c and large 'tHooft coupling 
expansions. In the holographic approach, infinite tower 



of color-singlet modes with mesonic quantum numbers as 
p, fli, p' , a[, /o", • • • , appears in the mode expansions of 
the five-dimensional gauge field Ap{x,z) on the D8 brane 
as in Eq.(7). In this sense, the baryon as the instanton 
on the five-dimensional gauge theory before the mode 
expansion is in principle to be composed by the infinite 
tower of such color-singlet modes with mesonic quantum 
numbers. Therefore, the instability of the instanton only 
with the DBI sector is often regarded that the infinite 
number of (axial) vector mesons p, oi, p', p" , 
would affect the low-energy soliton feature to give the 
zero size. Then, in order to describe the baryon as an in- 
stanton with finite size in the holographic approach, the 
inclusion of the Chern-Simons (CS) sector of the effective 
action of the D8 brane is claimed, whereas the CS sec- 
tor corresponds to the higher order contribution of the 
'tHooft coupling expansion relative to the DBI sector. 

However, as emphatically noted in Sec. II, the appear- 
ance of certain cutoff scale like A/kk ^ IGeV in holo- 
graphic approach should be essential to be dual of QCD. 
In fact, the mesonic mass spectra predicted from holo- 
graphic QCD starts to deviate from the experimental 
data beyond the Mkk scale, indicating that the holo- 
graphic duality with QCD is mainly maintained almost 
below Mkk- This tendency of the mass spectra also de- 
notes that such color-singlet modes beyond cutoff scale 
Mkk do not directly correspond to physical mesons in 
QCD. Therefore, one cannot manifestly conclude that 
the instability of the instanton only with the DBI sec- 
tor really arises from the physical effect in QCD. 

Furthermore, there also exist the infinite number of 
non-QCD modes, e.g., the Kaluza-Klein modes with large 
mass ^ 0{Mkk) through the Kaluza-Klein compactifi- 
cation of the D4 brane [3]. Therefore, if the baryon as 
the instanton is to be really composed by the infinite 
tower of the color-singlet non-QCD modes even beyond 
Mkk, there is no reason to cast away the Kaluza-Klein 
modes in the baryon analysis. In fact, such Kaluza-Klein 
modes still have the possibility to affect baryon proper- 
ties, giving, e.g., the stability of the instanton even with- 
out the CS sector of the effective action of the D8 brane, 
through it is not to be manifestly proved yet. As a whole, 
there still seems to remain puzzling conflicts between the 
truncated resonance model and instanton models for the 
baryon analyses in holographic QCD. 

By looking back the long history of hadron physics, 
baryons seem to have inherent difficulties to be described 
relative to mesons. For example, the quark model suc- 
ceeded in the systematic classification of the hadrons in 
terms of their valence quarks [65, 66]. However, while the 
dynamics of mesons as the two-body composites can be 
relatively easily analyzed by solving the Bethe-Salpeter 
equation [67], we have to treat the complicated Faddeev 
equation [68] to describe the dynamics of the baryons 
as the three-body composites. Same difficulties are also 
found in the lattice QCD studies, for instance, with re- 
spect to the inter-quark potential of mesons [69] and 
baryons [70]. In fact, the quark-antiquark potential as 
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the internal nature of the meson could be successfully 
measured m lattice QCD in 1980 [71], while it takes al- 
most twenty years after that to find the well measurement 
of three-quark potential with the appearance of gluonic 
Y-type flux tube as the internal nature of the baryon [72] . 
The difficulty of the baryons gets more apparent in the 
large- TVc QCD [8], which has provided a powerful pcrtur- 
bative treatment with the l/Nc expansion for the non- 
perturbative aspects of the strong interaction. In fact, 
large-iVc QCD achieved large successes in the explana- 
tions for the hadron phenomenology, e.g., the Okubo- 
Zweig-Iizuka (OZI) rule [50], the A/ = 1/2 rule [73], 
and the narrowness of meson resonances [8, 50]. It also 
provided a quantitative formula for the large r]' meson 
mass as the Veneziano-Witten formula [74] . On the other 
hand, the baryon docs not even appear as the dynami- 
cal degrees of freedom in the large- A^c QCD. because the 
baryon mass is proportional to 0{Nc). In other words, Nc 
quarks are needed to construct the SU(iVc) color-singlet 
composite as a baryon only from quarks, while a me- 
son can still be constructed as the two-body color-singlct 
composite from a quark and an antiquark. In the end, 
the baryon appears as a "non-local object", i.e., the soli- 
ton of the meson fields in the large- A^c QCD. As a whole, 
the baryons tend to suffer from the many-body difficul- 
ties relative to the mesons, by seeing the long hadron 
history. 

Now, we suppose here that the recent baryon analyses 
in holographic QCD may suffer from the same kind of 
difficulties inherent in the baryon itself. In fact, meson 
properties are successfully described in the framework 



of holographic QCD as the Sakai-Sugimoto model [3] 
whereas the baryons seem to be not sufficiently described 
yet. One should note that the baryon mass splitting cor- 
responds to the higher-order contributions of the large- 
Nc expansion. In terms of the superstring theory, such 
baryon mass splitting corresponds to the string loop ef- 
fect beyond the classical supergravity, which is fairly un- 
tractable. These naive considerations imply that the 
baryon still drags the essential difficulties even by start- 
ing from the superstring theory. We have to sincerely 
reconsider the meaning of the difficulty existing in the 
baryon itself and identify an origin of the problems, which 
would inspire us to cure the recent conflicts of the baryon 
analyses in the holographic approach for the future. 
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